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Abstract 

We describe the motivation behind the recent formulation of a nonperturbative 
path integral for Lorentzian quantum gravity defined through Causal Dynamical Tri- 
angulations (CDT). In the case of two dimensions the model is analytically solvable, 
leading to a genuine continuum theory of quantum gravity whose ground state de- 
scribes a two-dimensional "universe" completely governed by quantum fluctuations. 
One observes that two-dimensional Lorentzian and Euclidean quantum gravity are 
distinct. In the second part of the review we address the question of how to incorpo- 
rate a sum over space-time topologies in the gravitational path integral. It is shown 
that, provided suitable causality restrictions are imposed on the path integral histo- 
ries, there exists a well-defined nonperturbative gravitational path integral including 
an explicit sum over topologies in the setting of CDT. A complete analytical solu- 
tion of the quantum continuum dynamics is obtained uniquely by means of a double 
scaling limit. We show that in the continuum limit there is a finite density of infinites- 
imal wormholes. Remarkably, the presence of wormholes leads to a decrease in the 
effective cosmological constant, reminiscent of the suppression mechanism considered 
by Coleman and others in the context of a Euclidean path integral formulation of 
four-dimensional quantum gravity in the continuum. In the last part of the review 
universality and certain generalizations of the original model are discussed, providing 
additional evidence that CDT define a genuine continuum theory of two-dimensional 
Lorentzian quantum gravity. 
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1 Introduction 

Finding a consistent theory of quantum gravity which gives a fundamental quantum de- 
scription of space-time geometry and whose classical limit is general relativity lies at the 
root of our complete understanding of nature. However, nearly one century has passed 
by since Einstein's invention of general relativity, and still very little is known about the 
ultimate structure of space-time at very small scales such as the Planck length. Because 
of the enormous energy fluctuation predicted by the Heisenberg uncertainty relation, the 
geometry at such short scales will have a non-trivial microstructure governed by quantum 
laws. The difficulty of finding a description of this microstructure is on the one hand due 
to the lack of experimental tests and further complicated by the fact that four-dimensional 
quantum gravity is perturbatively non-renormalizable. This also holds for supergravity or 
perturbative expansions in the string coupling in string-theoretic approaches. One road to 
take is therefore to try to define quantum gravity nonperturbatively. 

Within those nonperturbative approaches of quantum gravity there are several attempts 
which suggest that the ultraviolet divergences can be resolved by the existence of a minimal 
length scale, commonly expressed in terms of the characteristic Planck length. A famous 
example is loop quantum gravity [TJ El El ; in this canonical quantization program the dis- 
crete spectra of area and volume operators are interpreted as an evidence for fundamental 
discreteness. Other approaches, such as four dimensional spin-foam models [I] or causal set 
theory jSlE], postulate fundamental discreteness from the outset. Unfortunately, neither 
of these quantization programs has succeeded so far in recovering the right classical limit. 

A recent alternative is Lorentzian quantum gravity defined through Causal Dynamical 
Triangulations (CDT). In this path integral formulation a theory of quantum gravity is 
obtained as a continuum limit of a superposition of space-time geometries. Thereby, in 
analogy to standard path integral formulations of quantum mechanics and quantum field 
theory, one uses an intermediate regularization through piecewise flat simplicial geome- 
tries. 1 

Causal dynamical triangulations were first introduced in two dimensions [7J |H] as a 
nonperturbative path integral over space-time geometries with fixed topology which can 
be explicitly solved, leading to a continuum theory of two-dimensional Lorentzian quantum 
gravity. Later it was extended to three CHI HI] and four dimensions (01 H21 UH1 UH HHj 
(see also ^Bj for a general overview). Among those recent developments there are very 
promising results regarding the large scale structure of space-time fF2\ ITS*] : Firstly, the 
scaling behavior as a function of space-time volume is that of a genuine isotropic and 
homogeneous four-dimensional universe; this is the first step towards recovering a sensible 
classical limit. Moreover, after integrating out all dynamical variables apart from the scale 
factor a(r) as a function of proper time r, it describes the simplest minisuperspace model 
used in quantum cosmology. 

Very interesting is the analysis of the microstructure of space-time: recent numerical 



however, whereas in quantum theory the piecewise straight paths are imbedded in a higher-dimensional 
space, this is not the case for the piecewise flat geometries whose geometric properties are intrinsic. 
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results showed [HI HH] that in this setting space-time does not exhibit fundamental dis- 
creteness. Instead, one observes a dynamical reduction of the dimension from four at large 
scales to two at small scales. This gives an indication that nonperturbative Lorentzian 
quantum gravity provides an effective ultraviolet cut-off through a dynamical dimensional 
reduction of space-time. 

In this review we give a pedagogical introduction into CDT in two dimensions. The sim- 
ple structure of two-dimensional gravity serves as a good playground to address fundamen- 
tal concepts of the model which are also of significant relevance for the higher- dimensional 
realizations of CDT. 

One interesting fundamental question we want to address here is whether a sum over 
different space-time topologies should be included in the gravitational path integral. Since 
topology changes [T7] naturally violate causality, the sum over topologies is usually con- 
sidered in Euclidean quantum gravity, where this issue does not arise. However, even in 
the simplest case of two-dimensional geometries the number of configurations contributing 
to the path integral grows faster than exponentially which makes the path integral badly 
divergent. Various attempts to solve this problem in the Euclidean context have so far 
been unsuccessful [TH], which leaves this approach at a very unsatisfactory stage. 

Other attempts to define a path integral including a sum over topologies are either 
semiclassical or assume that a certain handpicked class of configurations dominates the path 
integral, without being able to check that they are saddle points of a full, nonperturbative 
description ^H]- However, we are not aware of any nonperturbative evidence supporting 
these ideas. 2 

In the context of two-dimensional CDT it has recently been shown that one can unam- 
biguously define a nonperturbative gravitational path integral over geometries and topolo- 
gies (2HI22] 3 - The idea of how to tame the divergences is to use certain causal constraints 
to restrict the class of contributing topologies to a physically motivated subclass. In the 
concrete setting these are geometrically distinguished space-times including an arbitrary 
number of infinitesimal "wormholes", which violate causality only relatively mildly. By 
this we mean that they do not necessarily exhibit macroscopic causality violations in the 
continuum limit. According to observable data this is a physically motivated restriction. 
This makes the nonperturbative path integral well defined. As an interesting result one 
observes that in the resulting continuum theory of quantum gravity the presence of worm- 
holes leads to a decrease in the effective cosmological constant. This connects nicely to 
former attempts to devise a mechanism, the so-called Coleman's mechanism, to explain the 
smallness of the cosmological constant in the Euclidean path integral formulation of four- 
dimensional quantum gravity in the continuum in the presence of infinitesimal wormholes 

Another interesting aspect is the discussion of universality of two-dimensional 
Lorentzian quantum gravity defined through CDT. We provide explicit evidence for uni- 

2 For another attempt to include the sum over topologies in the context of three-dimensional quantum 
gravity the reader is referred to |2(J| . 

3 Also see [231 |^j for additional information. 
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versality of the model which ensures that CDT defines a genuine continuum theory of 
two-dimensional quantum gravity and does not depend on the details of the regularization 
procedure. 

The remainder of this review is structured a follows. In Section El we give an intro- 
duction to two-dimensional Lorentzian quantum gravity defined through CDT. We show 
that it gives rise to a well-defined continuum theory of two-dimensional quantum grav- 
ity, whose time evolution is unitary. Further, the restriction to the causal geometries in 
the path integral leads to a physically interesting ground state of a quantum "universe" 
which also shows that quantum gravity with Euclidean and Lorentzian signature are dis- 
tinct theories. In Section |3] we address the question of how to incorporate the sum over 
topologies in the gravitational path integral. After restricting the contributing topologies 
to a physically motivated subclass including arbitrary numbers of infinitesimal wormholes, 
the path integral is well-behaved. The continuum theory can be obtained by means of a 
well-defined double scaling limit of the couplings. We observe a finite density of worm- 
holes which leads to a decrease in the effective cosmological constant. In Section 0] we 
show evidence for universality of the two-dimensional model. Further, we develop several 
one-to-one correspondences between two-dimensional CDTs and certain one-dimensional 
combinatorial structures, such as heaps of dimers and random walks. Appendices lATIDl 
provide supplementary information for Sectional In Appendix El we give a brief summary 
of the results on Lorentzian angles which naturally appear in the framework of CDT. In 
Appendix IB1 triangulations with different boundary conditions are discussed. In Appendix 
lUl we give an alternative derivation of the continuum quantum Hamiltonian obtained in 
Section |21 Regarding the calculation of continuum quantities further details are given in 
Appendix \D\ In Appendix |E] we discuss other double scaling limits related to the model 
in Sectional which we discarded as unphysical. 



2 2D Lorentzian quantum gravity 

In this section we give an introduction to two-dimensional Lorentzian quantum gravity de- 
fined through causal dynamical triangulations (CDT). The section is structured as follows: 
In Section 12.11 we introduce the reader to the concepts and problems of how to define a 
nonperturbative path integral for quantum gravity. In Section 12.21 an explicit definition 
of the gravitational path integral as a regularized sum over space-time triangulations is 
presented. The problem of how to incorporate Lorentzian space-times in this framework 
and how to give a well-defined prescription for a Wick rotation is discussed in Section 12.31 
this leads to the notion of causal dynamical triangulations. In Section |2~H a complete ana- 
lytic solution to the discrete problem is presented. Performing the continuum limit of the 
discretized model, as done in Section l2~ol leads to a continuum theory of two-dimensional 
Lorentzian quantum gravity, which we will interpret in terms of its physical observables in 
Section I2~fi1 
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2.1 A nonperturbative path integral for quantum gravity 

In this section a notion of a nonperturbative path integral for quantum gravity is intro- 
duced. In contrast to usual path integral formulations in quantum mechanics and quantum 
field theory the gravitational path integral is more involved due to the diffeomorphism- 
invariance of the theory, which in terms of local coordinate charts are smooth invertible 
coordinate transformations x^ \— > y fl (x tl ), and due to the absence of a preferred metric 
background structure. 

The aim of this section is to expose this problem and to give first indications on how 
to define a gravitational path integral. Most of the arguments presented in this section 
can be found in |2Zj . For a more comprehensive account of path integrals in quantum field 
theory and the connection to critical phenomena the reader is referred to |28j . 

Before going into the conceptual details of the problems one encounters when defin- 
ing a gravitational path integral, let us first recall the path integral representation of an 
one-dimensional (non-relativistic) quantum mechanical problem described by the time de- 
pendent Schrodinger equation 

-i^|<M> = £(t)IM- (i) 

The time evolution of the wave function is given by 

|V,t"> = w(tV) \ip,t'), U(t,t) = l, (2) 



where U(t,t') is defined by H(t) = -i d t U(t,t) 
time evolution, (J2J), can be written as 



. In the position representation, the 



ijj(x", t") = / dx' G(x", x'; t", t') ^(x', £'), (3) 



where G(x",x';t",t') = (x" \ U(t",t r ) \ x') is the so-called propagator or Feynman kernel. 
For the case that the Hamiltonian H(t) is bounded from below, U(t,t') satisfies the semi- 
group property 

U{t", t') = U{f, t) U{t, t'), t' ^t^ t". (4) 

This property (jlj allows us to define the propagator as a limiting procedure of products 
of evolution operators corresponding to infinitesimal time intervals e = (t" — t')/N (Figure 



N-l . N-l 

G(x" ,x';t" ,t') = limA~ N JJ / dx k JJ <x i+i | U(t j+1 ,tj) \xj) 

k=l J j=0 
N-l „ ( N-l 



lim A N \\ J dx k exp < i ^ e£(x j+1 , (x j+ i - Xj)/e) > , (5) 



k=l " (. j=0 
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Figure 1: Illustration of the path integral for a one-dimensional non-relativistic quantum 
mechanical problem, e.g. a propagating particle. One possible path of the configuration 
space (path space) is drawn. The "virtual" particle is propagating from x to ijv in a 
piecewise linear path of N steps of time e = (£" — t')/N each. 



where A is a normalization factor and £ is the Lagrangian corresponding to H(t). The 
integral in (jSJ) is taken over piecewise linear paths x(t) from x' = Xq to x" = x^. The right 
hand side of (jSJ) is called the path integral and is often written symbolically as 

G(x",x';t",t') = [ V[x(t)}e zS[x ^\ (6) 
Jv 

where is a functional measure on the path space V [2H|- The weight of each path 

is given by the classical action S[x(t)}. 

Note that © is mathematically not well-defined. To actually perform the integral in 
(JHJ) one usually does an analytic continuation on the time variable t from the real to the 
imaginary axis, i.e. 1 1— ► r = it. This so-called Wick-rotation makes the integrals real (201 
by mapping i S[x(t)] i— > —S eu [x(t)], where S eu [x(t)} is the Euclidean action. After having 
performed the integration in the Euclidean sector one applies the inverse Wick rotation 
t i— > —it to get back to the physically meaningful results. 

The question which arises now is: Can one define a path integral for quantum gravity 
using a similar strategy like the one used to arrive at (JBJ)? The first step to answer this ques- 
tion is to find the analogue of path and the path space in a gravitational sense. One might 
think that the equivalent of the path space V in the case of the one- dimensional quantum 
mechanical problem is the space of all geometries in the gravitational case. To clarify this, 
let us become more precise and provide a definition of these concepts. Classically, a geom- 
etry is a space-time (M,g), a smooth manifold M equipped with a metric tensor field 
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with Lorentzian signature [HI]. Due to the diffeomorphism invariance of the theory any 
two metrics are equivalent if they can be mapped onto each other by a coordinate trans- 
formation. Therefore all physical degrees of freedom are encoded in the equivalence class 
\9tw] = 9[iv/Di&(M) which we simply refer to as a geometry. From this follows a natural 
notion of the space of all geometries as the coset space Geom(M) = Metrics (M)/Diff(M). 




Figure 2: Illustration of the gravitational path integral; one geometry of the configuration 
space is drawn. The geometry [g^ u ] has a time-sliced structure with respect to proper time 
£; it has an initial spatial geometry [g'^] at proper time t' and final spatial geometry [g 1 ^] 
at proper time t". Note that in the case of two dimensions the spatial hypersurface E is 
one-dimensional and the spatial geometry [gry] is totally characterized by its length L. 

Having defined the space of all geometries we can use a foliation of the space of geome- 
tries to schematically define the gravitational path integral from an initial spatial geometry 
[g'ij] at proper time t' to a final spatial geometry [g'/j] at proper time t" (Figure |2J) by 

GW].W^.O= / v[g»Ae iSM - (7) 

JGeom(M) 

Here, T>[g^ u ] is the functional measure on the space of all geometries Geom(M) and each 
geometry is weighted by the classical Einstein-Hilbert action Seh[sVm']- 4 Before trying to 
give a concrete meaning to (JTj) in the next section, let us first discuss the conceptional 
problems of how to define and actually evaluate (JJJ) . 

First of all, T>[g^ u ] has to be defined in a covariant way to preserve the diffeomorphism 
invariance. Since there is no obvious way to parametrize geometries, one would have to 

4 We assume that each geometry is weighted by the classical Einstein-Hilbert action including a term 
with a positive cosmological constant. One might also include higher curvature terms in the action. 
However, as we will see in Section 14.11 the resulting continuum theory belongs to the same universality 
class as the continuum theory obtained from Seh- 



2.2 Geometry from simplices 
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introduce covariant metric field tensor. Unfortunately, to perform further calculations one 
would then have to gauge fix the field tensors which would give rise to Faddeev-Popov 
determinants j22] whose nonperturbative evaluation is exceedingly difficult. 5 

The second problem is due to the complex nature of the integrand. In quantum field 
theory this problem is solved by doing the Wick rotation t i— > r = it, where t is the time 
coordinate in Minkowski space (t, x). Clearly, a prescription like this does not work out in 
the gravitational setting, since all components of the metric field tensor depend on time. 
Further t i— > r = it is certainly not diffeomorphism invariant (as a simple example consider 
the coordinate transformation t i— > t 2 for t > 0). Hence, the question arises: What is the 
natural generalization of the Wick rotation in the gravitational setting? 

Finally, since we are working in a field theoretical context, some kind of regularization 
and renormalization will be necessary, and again, has to be formulated in a covariant way. 
In quantum field theory the use of lattice methods provides a powerful tool to perform 
nonperturbative calculations, where the lattice spacing a serves as a cutoff of the theory. An 
important question to ask at this point is whether or not the theory becomes independent 
of the cutoff. Consider for example QCD and QED on the lattice. All evidence suggests 
that QCD in four space-time dimensions is a genuinely continuum quantum field theory. 
By genuine continuum quantum field theory we mean a theory which needs a cutoff at 
an intermediate step, but whose continuum observables will be independent of the cutoff 
at arbitrary small scales. In QED the situation seems to be different, one is not able 
to define a non-trivial theory with the cut-off removed, unless the renormalized coupling 
e-ren = (trivial QED). Therefore, QED is considered as a low-energy effective theory of 
an elaborate theory at the Planck scale. However, as we will in Section 12 .5[ for the case 
of the gravitational path integral the continuum limit exists and the resulting theory is 
independent of the cutoff, in the same way as QCD is. 

2.2 Geometry from simplices 

As motivated in the previous section, regularization through lattice methods might be a 
powerful tool for a nonperturbative path integral formulation of quantum gravity. This 
obviously requires some kind of discretization of the geometries. 

The use of discrete approaches to quantum gravity has a long history [3^1 • In classical 
general relativity, the idea of approximating a space-time manifold by a triangulation of 
space-time goes back to the early work of Regge j2Zj and was first used in |3H1 EH] to 
give a path integral formulation of gravity. By a triangulation we mean a piecewise linear 
space-time obtained by a gluing of simplicial building blocks. This one might think of as 
the natural analogue of the piecewise linear path we used to describe the path integral 
of the one-dimensional quantum mechanical problem. In two dimensions these simplicial 
building blocks are flat Euclidean or Minkowskian triangles, where flat means isomorphic 

5 See |33l I34j for an evaluation in the setting of two-dimensional Euclidean quantum gravity in the 
light-cone gauge. In |35| a calculation for three- and four-dimensional Lorentzian quantum gravity in the 
proper-time gauge is presented. It is anticipated there that the Faddeev-Popov determinants cancel the 
divergences coming from the conformal modes of the metric nonperturbatively. 
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2 2D LORENTZIAN QUANTUM GRAVITY 



to a piece of Euclidean or Minkowskian space respectively. One could in principle assign a 
coordinate system to each triangle to recover the metric space (M,g^ u ), but the strength 
of this ansatz lies just in the fact that even without the use of coordinates, each geometry 
is completely described by the set of edges length squared {If} of the simplicial building 
blocks. This provides us with a regularized parametrization of the space of all geometries 
Geom(M) in a diffeomorphism invariant way, and hence, is the first step towards defining 
the gravitational path integral (0). 

For the further discussion it is essential to understand how a geometry is encoded in 
the set of edges length squared {If} of the simplicial building blocks of the corresponding 
triangulation. In the case of two dimensions this can be easily visualized, since the trian- 
gulation consists just of triangles. Further, the Riemann scalar curvature R(x) coincides 
with the Gaussian curvature K(x) up to a factor of 1/2. There are several ways to reveal 
curvature of a simplicial geometry. The most convenient method is by parallel transporting 
a vector around a closed loop. Since all building blocks are flat, a vector parallel trans- 
ported around a vertex v always comes back to its original orientation unless the angles 8{ 
of the surrounding triangles do not add up to 2n, but differ by the so-called deficit angle 
e v = 27r — J2i^ v @* (Figure BJ). The Gaussian curvature located at the vertex v is then given 
by 

K v = (8) 

' V 

where V v is the volume associated to the vertex v; more precisely, the dual volume of the 
vertex. Note that the curvature at each vertex takes the form of a conical singularity. 
Then one can write the simplicial discretization of the usual curvature and volume terms 
appearing in the two-dimensional Einstein-Hilbert action, 

- / d 2 xy/\detg\R(x) -> Ve„, (9) 
I d 2 aV|det0| -> J^K, (10) 

where 71 = {If} is a triangulation of the manifold M described by the set of edge length 
squared. From this one can write down the simplicial discretization of the two-dimensional 
Einstein-Hilbert action, the so-called Regge action, 



where k is the inverse Newton's constant and A the cosmological constant. It is then an 
easy exercise of trigonometry to evaluate the right hand side of (jllj) in terms of the squared 
edge lengths. 

The treatment explained above can also be generalized to an arbitrary dimension d in 
a straightforward way [40J. 

From this notion of triangulations one can now write the gravitational path integral as 
the integral over all possible edge lengths j 7S>1, where each configuration TZ is weighted by 



2.3 Lorentzian nature of the path integral 
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Figure 3: Illustration of a positive (a) and negative (b) space-like deficit angle e v at a 
vertex v. 



the corresponding Regge action (jllj) . A potential problem one encounters with this ansatz 
is an overcounting of possible triangulations, due to the fact that one can continuously 
vary each edge length. Further, one still has to introduce a suitable cut-off for the length 
variable I. Among other things, this motivated to the approach of "rigid" Regge calculus or 
dynamical triangulations where one considers a certain class T of simplicial space-times as 
an explicit, regularized version of Geom(M), where each triangulation TgT only consists 
of simplicial building blocks whose space-like edges have all the same edge length squared 
I 2 = a 2 and whose time-like edges have all the same edge length squared I 2 = —a 2 . Here, 
the geodesic distance a serves as the short-distance cutoff, which will be sent to zero later. 
At this point it is important to notice that fixing the edge length squared is not a restriction 
on the metric degrees of freedom. One can still achieve all kinds of deficit angles of either 
sign or directions by a suitable gluing of the simplicial building blocks. 

From this one is able to give a definite meaning to the formal continuum path integral 
((2j) as a discrete sum over inequivalent triangulations, 



where 1/C(T) is the measure on the space of discrete geometries, with C{T) = |Aut(T)| 
the dimension of the automorphism group of the triangulation T. 

2.3 Lorentzian nature of the path integral 

After we have given an explicit meaning to gravitational path integral as a sum over 
inequivalent triangulations ()12j). there still remains the problem of how to actually perform 
the sum, since it is up to now unclear what the prescription of a Wick rotation looks like, 
i.e. how we are going to analytically continue to the Euclidean sector of the theory. 




(12) 
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To avoid this problem one first considered dynamical triangulations jHJEj on Euclidean 
geometries made up of Euclidean simplicial building blocks, where one does the ad hoc 
substitution 



where one simply replaced by hand the Lorentzian geometries with their Euclidean coun- 
terparts. The reason for doing so was not necessarily that one thought that Euclidean 
manifolds are more fundamental than Lorentzian manifolds, it was just that one had no a 
priori prescription for a Wick rotation. 

The potential problem with the substitution (JT3Jl is that one integrates over geometries 
which know nothing about time, light cones and causality. Hence, there is no a priori 
prescription of how to recover causality in the full quantum theory. It is unlikely that this 
can be done by just performing the inverse Wick rotation t i— > —it. Another thing which 
fails is that Geom(M e ") has too many configurations and one integrates over a large class 
of highly degenerate acausal geometries. This problem exists in higher dimensions d ^ 3, 
where it causes the absence of a well defined continuum theory. More precisely, due to the 
degenerate geometries, Euclidean dynamical triangulations in d = 3, 4 has an "unphysical" 
phase structure, where neither of the phases seems to have a ground state that resembles 
an extended geometry. 6 

This led to the approach of Lorentzian dynamical triangulations or causal dynamical 
triangulations (CDT) which was first introduced in two dimensions in [7j, further elabo- 
rated in [HI El and has recently produced interesting results in d = 4 [T5J CH H21 EE] as 
presented in the introduction. The main point of CDT is to insist in taking the Lorentzian 
structure, i.e. the inherent light-cone structure, seriously from the outset. Due to the 
triangulated structure one is able to find a well-defined Wick rotation on the full path 
integral (|12j) as we will discuss in the following. Once having defined the Wick rotation all 
calculations, including the continuum limit, are performed in the Euclidean sector, where 
at the end one goes back to the physical signature by performing the inverse Wick rotation 
(in the continuum theory). 

Due to the inherent causal structure, the Lorentzian model is generally inequivalent 
to the Euclidean model. 7 One might think of the causal structure as being equivalent 
to choosing a measure on the path integral (fT2Jl which suppresses acausal geometries and 
therefore also some of the highly degenerate geometries. 

To be able to define the Wick rotation in this context, let us first clarify which are the 
causal triangulations T that contribute to the sum in the path integral (|12j). 

We consider "globally hyperbolic" simplicial manifolds with a sliced structure, where 
the one-dimensional spatial hypersurfaces have fixed topology of S 1 ; this will lead to a 
de-Sitter space-time in the continuum theory. Further we will not allow for topology 
changes of the spatial slices; in how far topology changes can be incorporated within this 

6 See I3U and references therein for a detailed discussion of the phase structure of four-dimensional 
dynamical triangulations. 

7 In two dimensions the Euclidean model can be related to the Lorentzian one by a non-trivial procedure 
of integrating out baby universes [5]. 




(13) 
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Figure 4: Illustration of the causal structure in a flat triangulated manifold. The triangles 
are glued according to the sliced structure. Due to the inherent causal structure of the 
Minkowskian triangles, the whole simplicial manifold has a causal structure. 



framework will be discussed in Section El The sliced structure of the simplicial manifold 
provides us with a preferred notion of "proper time", namely, the parameter labeling 
successive spatial slices. Note that this use of time is not a gauge choice, since proper 
time is naturally defined in a diffeomorphism invariant way. Several spatial slices are then 
connected by Minkowskian triangles 8 with one space-like edge of length squared I 2 = +a 2 
and two time-like edges of length squared I 2 = —a 2 . Since all simplicial building blocks are 
Minkowskian triangles which have an intrinsic light-cone structure, we have local causality 
relations which, if one uses certain gluing rules, imply a global causal structure on the 
whole manifold (Figure EJ). 

From this notion, one is now able to define the Wick rotation on a very elementary 
geometric level. By virtue of the triangulated structure one can define a Wick rotation on 
each simplicial building block respectively by performing a Wick rotation on all time-like 
edges It of the Lorentzian triangles, i.e. It i— > il t . The Wick rotation acting on the whole 
simplicial manifold T is then defined by the following injective map 

W : T lor = {T, I 2 = a 2 , l 2 t = -a 2 } h-> T eu = {T, l 2 s = a 2 , 1 2 = a 2 }. (14) 

In Appendix |X] it is shown how to implement W in a mathematically clean way. Further, 
it is shown there that in the path integral, the Wick rotation W implements precisely the 
desired analytic continuation on the Regge action, 

W . giWeC^) ^ g-SawM. (15) 

Now we have all ingredients in hand to define the path integral of the regularized 
model. As we will see in the following sections, evaluation of the propagator then reduces 
to a simple combinatorial problem and applications of the theory of critical phenomena. 

8 A brief description of how to assign angles and volumes to Minkowskian triangles is given in Appendix 

El 
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t + 2 




t + l 



t 



9 

+a 



Figure 5: Section of a 2d Lorentzian triangulation consisting of space-time strips of height 
At = 1. Each spatial slice is periodically identified, such that the simplicial manifold 
has topology [0, 1] xS 1 . One sees that a single strip with lower boundary length l\ and 
upper boundary length l 2 consists exactly of h up pointing triangles and l 2 down pointing 
triangles. 



2.4 Discrete solution: The transfer matrix 

In this section we evaluate the gravitational path integral of the two dimensional CDT 
model following the procedure explained above [Zj. 

The space-time manifolds we are considering consist of just one type of simplicial build- 
ing blocks, namely, flat Minkowskian triangles with one space-like edge of length squared 
l 2 s =a 2 and two time- like edges of length squared I 2 = —a 2 . Each triangle then has a volume 
of ^a 2 (Appendix [X]). Our globally hyperbolic simplicial manifold is built up of space-time 
strips of height At = 1, where each strip consists of up and down pointing Minkowskian 
triangles (Figure |3J). The foliation parameter t is interpreted as the discretized version of 
"proper" time T = a ■ t. Each spatial slice at time t is chosen to have periodic boundary 
conditions with fixed topology S 1 of the sphere. (Different boundary conditions are dis- 
cussed in Appendix IE1IU|) . Note that each spatial geometry is completely characterized by 
its discrete length /, where L = a ■ I is the spatial length. Further, due to the S 1 topology, 
the symmetry factor of a spatial slice of length I is just C(l) =1. 

After we have defined the triangulations contributing to the path integral, the next 
step is to determine their corresponding weights, i.e. the Regge action. Recall that the 
Einstein-Hilbert action in two dimensions is given by 



where K = G N is the inverse Newton's constant. A simplification occurs due to the 
topological character of the curvature term in two dimension, 




(16) 




(17) 
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where x(M) = 2 — 2g is the Euler characteristic of the manifold M and q is the genus of 
M. Hence, since we are not allowing for spatial topology changes, the exponential of % 
times this term is just a constant phase factor which does not influence the dynamics and 
therefore can be pulled out of the path integral. If one actually allows for topology changes 
this term will be more involved as we will see in Section El Since all triangles have the 
same volume, the Regge action fTTf takes the following simple form 



S Reggc (T) = Xa 2 N(T), (18) 

where A is the bare cosmological constant and N(T) the number of triangles in the tri- 
angulation T. Note that a factor of ^ coming from the volume term has been absorbed 
into A. Using fT2|) and (|18jl . the discrete gravitational path integral over the set of two 
dimensional causal simplicial manifolds T with an initial boundary U n and a final boundary 
l out , consisting of t = t out — t in time steps, can be written as 

Gf{l m Xuut)= ±e**»M (19) 



TeT: 



For simplicity we can remove the symmetry factor \jli n by marking a vertex on the initial 
boundary, hence define 

G\(lini lout: t) hnG^ (lini louti t) ■ (20) 

The unmarked propagator can be easily recovered at a later stage. 

The next step of the calculation will be to Wick rotate G x (l in , l out ; t) to the Euclidean 
sector. Following the prescription for the Wick rotation W as defined in (|14j) . we get 

G x (l m Xuut)= J2 e*~ Xa2N ^^G x (l m Xuut)= Yl e-' a2 ^ T \ (21) 

TeT : TeW(T): 

kn^lout k„—*lout 

where A and A differ by a 0(1) constant, due to the volume difference of the Euclidean 
and Minkowskian triangles (Appendix EJ. 

An essential quantity for the solution of the model is the transfer matrix T which 
contains all dynamical information of the system. The transfer matrix can be represented 
by its matrix elements with respect to the initial state | li n ) and the final state | l ou t), i.e. 
the kernel 

G \{lini lout: t 1) (lout \ T | ij n ) , (22) 

which is often referred to as the one-step propagator. From the semi-group property of the 
propagator, 

G x (h,l2',tx + t 2 ) = ^GxihJ^Gx&h-M, (23) 

i 

Gx{h,WM + l) = Y, G ^,l-M)Gx(l,WA), (24) 
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it follows that the propagator G\(li n , l ou u t) for an arbitrary time t is obtained by iterating 
(I22[) t times, yielding 

G\(lin,lout',t) = (li n \ T* | l ut)- (25) 

In the continuum language this expression is closely related to the continuum partition 
function, 

Z T = lim Trf*, (26) 

a— >0 

where in the limit the product T = t-a is kept fixed. One sees that knowing the eigenvalues 
of the transfer matrix is the key to solving the general problem. 

Let us therefore concentrate on triangulations T of just one single time step: For such a 
strip with initial boundary of length U n and final boundary of length l out the total number 
of triangles is just N(T) = l in +l out (cf. FigureEJ). Evaluating the eigenvalues of the transfer 
matrix then reduces to a simple counting problem, 

G x (l m J out -t = l) = e- Xa2 ^ +l °^ Yl l = e^ a2 ^ +l ^-^—( hn ^ lout ). (27) 

TeW(T)- * ^ ' 

For later purposes it is convenient to introduce the generating function for the propa- 
gator Gx(hn,lout',t), defined as 

G(x,y;g;t) = ^ G x {l in ,l out ;t)x l ^y 1 ^, (28) 

where we used the notation g = e~ Xa2 . This expression can also be used to rewrite the 
semi-group property or composition law (J24)) . yielding 

/dz 
——G(x,z~ 1 ;g;t 1 )G(z,y;g;t 2 ). (29) 
Lniz 

The quantities x and y can be seen as purely technical tools of the generating function 
formalism, but one can also view them as boundary cosmological constants 

This interpretation is useful in the context of renormalization as we will see in the next 
section. The generating function of the one-step propagator can be obtained from ()27|) and 
(j2H| by usual techniques, yielding 9 

2 

G(x, y; g; 1) = (31) 

(1 - gx)(l - g{x + y)) 



9 Note that the asymmetry in x and y is due to the marking in the initial spatial boundary. "Unmarking" 
at a later stage will recover symmetry. 
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The joint region of convergence of (|<U|) is given by 

\x\ < 1, \y\ < 1, |#| < -. (32) 

Another way of obtaining fpTT^ is by graphical methods. Associating a factor of gx to 
every up pointing triangle "A" and a factor of gy to every down pointing triangle "v'\ 
the sum over all possible triangulations of one strip with circular boundary conditions can 
be written as 

G(x,y;g;l) = £J + + ^7 + ... 

oo / oo \ 00 

= E ^Eow) 1 -E(^) fc - ( 33 ) 

fc=0 \ «=0 / k=0 

The subtraction of the last summand has been performed to remove the degenerate cases 
where either entrance or exit loop have length zero. Evaluation of (|33|) shows the equiva- 
lence to (jSU>- 



2.5 Continuum limit 

In this section we want to construct the continuum limit of the discretized model, by 
making use of standard techniques of the theory of critical phenomena. Note that all the 
calculations are performed in the Euclidean sector of the theory. 

In terms of critical phenomena the continuum limit is obtained by fine-tuning the 
coupling constants, generally denoted by A, to the critical point A c of the phase transition. 
This means a divergence of the correlation length £ in the scaling limit 

£(A) ~ 77-^7 ^ °°, (34) 



A-A c| 



while the cutoff a goes to zero 



a(A) ~ |A- A c |" ^>0, (35) 

where z/>0 is a critical exponent. In terms of the cutoff a and the correlation length £ the 
physical correlation length is given by L = £(A) ■ a(A) and the continuum limit is defined as 
the simultaneous limit 

a(A)->0, £(A)^oo for fixed L = £(A) ■ a(A). (36) 

Having defined the general scheme to perform the continuum limit, let us now come back 
to the concrete problem of the discretized model described by the one-step propagator (J3TJ). 

From standard techniques of renormalization theory, we expect the couplings with 
positive mass dimension, i.e. the cosmological constant A and the boundary cosmological 
constants Aj n and X out to undergo an additive renormalization, 

A = ^ + A, A m = ^ + X, \ out = ^ + Y, (37) 
a z a a 
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where A, X and Y denote the corresponding renormalized values. Introducing the critical 
values 

g c = e~ Ax , Xc = e~ Ax in i y c = e ~ Ax out : (38) 

it follows that 

g = g c e- a2A , x = x c e~ aX , y = y c e~ aY . (39) 

The continuum limit can be performed by fine-tuning to the critical values 
(x,y,g) i — ► (x c ,y c ,g c ), as explained above. Further, it can be shown that the only 
sensible continuum theory can be obtained by choosing (x c , y c , g c ) = (1, 1, |) or (x c , y c , g c ) = 
(— 1, — 1, — |), which both lead to the same continuum theory. Hence, performing the con- 
tinuum limit can be done by simultaneously fine-tuning to the critical values with use of 
the following scaling relations: 

9 = \(l-a 2 A) + 0(a 3 ), (40) 

x = l-aX + 0(a 2 ), (41) 
y = 1-aY + 0(a 2 ). (42) 

Let us now show how we can reveal the continuum Hamiltonian in this scaling limit. 
First note that, as we have rewritten the composition law in terms of generating functions, 
i.e. (J29"j) . we can also rewrite the time evolution of the wave function (the analogue of (j3J)) 
in a similar manner, 

/dz 
G{x,z~ X ] l)if)(z, t). (43) 
2mz 

Introducing the scalings (|4L) p -()42 p and t=— into and expanding both sides to order a 
gives 



(l-aH{X 1 d x ) + 0{a 2 )) ip(X) 



io ° dZ f / 1 2X 2 -AXZ + Z 2 + 2K\ 9 A , s 

where we have defined ip{X) =ip(x = 1 — aX) and respectively for Z. For the expansion 
of the left hand side of (J44j) we used the formula for an infinitesimal time evolution of the 
wave function which is also used in the fundamental definition of the path integral, 

if?(X, T + a) = e- a " (x ' dx) ip{X, T). (45) 

Notice that the first term on the right hand side of (J43|) . ^3^, is the inverse Laplace 
transformed delta function 5(Z — X) as expected. Performing the integration in (J44j) . one 
can extract the Hamiltonian 



(46) 
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It is important to notice at this point that the Hamiltonian (|4fij) does not depend on 
possible higher order terms of the scaling relations (|4U |) -(|42 |) . One might check this by 
explicitly introducing higher order terms in the scaling relations, hence 

x = 1-aX + ^a 2 X 2 + 0{a 3 ), (47) 

y = 1-aY + ^a 2 Y 2 + 0{a 3 ). (48) 

Following the above procedure with the new scaling relations (j47| ) -()48 p . one sees that the 
resulting Hamiltonian does not depend on 7 and is still described by (|4"fij). 

Having obtained the effective quantum Hamiltonian in "momentum" -space X, one can 
go back to the Hamiltonian in length-space L by performing an inverse Laplace transfor- 
mation on the wave functions, 

/ io ° fix 
— e XL ^(X)- (49) 
■too 2m 

The inverse Laplace transformation is the continuum counterpart of the defining equation 
for the generating functions, formula (|28j). which one can think of as a discrete Laplace 
transformation. The effective quantum Hamiltonian for the marked propagator then reads 

H(L, ^) = -L— + 2AL, (marked). (50) 

This operator is selfadjoint on the Hilbert space 7i = £ 2 (M + , L~ 1 dL) . The effective quantum 
Hamiltonian corresponding to the unmarked propagator can then be obtained as 

dl) = ~ L M2 ~ 2 dL + 2AL ' ( unmarked )' ( 51 ) 

which is selfadjoint on the Hilbert space TC = £ 2 (R + , LalL). Notice that both Hamilton 
operators describe the same physical system. 10 In the following we want to analyze the 
effective quantum Hamiltonian (|51)1. calculate the spectrum, the eigenfunctions and from 
this further quantities like the partition function and the finite time propagator. A detailed 
calculation of these results is presented in Appendix [Dl 

The Hamiltonian consists of a kinetic term which depends on the spatial length of 
the "universe" and a potential term which depends on the (renormalized) cosmological 
constant. Further, it is important to notice that the Hamiltonian is bounded from below 
and therefore leads to a well-defined quantum theory. The spectrum of the Hamiltonian is 
discrete with equidistant eigenvalues, 



E n = 2V2A(n + 1) n = 0,1,2,.. (52) 



10 Consider the symmetric matrix element for the Hamiltonian of the marked case L^ 1 dLi(Li \ H \ L\). 
From (|60() it follows that unmarking can be done by absorbing a factor L^ 1 into the wave function, 
i.e. I L\) 1— >| Lf) = L^ 1 I L\). The unmarked Hamiltonian can then be obtained by evaluating 
dL^dLiiL^ H \Lx) = LidL x {Ll \ H u |L"). 
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The corresponding eigenfunctions are given by 

ip n (L) = Ae^ L iF^-n; 2; 2v / 2AX), dfi(L) = LdL, (53) 

where ra, fi + l;z) are the confluent hypergeometric functions. In this case the power 

series of the confluent hypergeometric function is truncated to a polynomial of degree n, 
namely, the generalized Laguerre polynomials, here denoted by L£(z), 

r(n + l)r(u + l) ! . . 

LJ (// + !)(// + 2)...(// + fc) ¥ " (55) 



fc=0 



From this one can see that the eigenfunctions ip n {L) form an orthonormal basis of the 
Hilbert space 7i, where the normalization factor is given by 

A n = 2^2A(n+ 1). (56) 

Having obtained the spectrum of the Hamiltonian one can easily calculate the (Euclidean) 
partition function by using 

oo 

Z T = Tr(e- TA ) = J2^ TEn - (57) 



n=0 



Hence, inserting (fo^j) into (|5T|) yields 



— 

e 2 v / 2AT 



mf^l = \ HMv^T) - l) . (5f 



Another interesting quantity to calculate is the finite time propagator or loop-loop corre- 
lator which is defined by 

G A (L 1} L 2 ;T) = (L 2 \e- TA \L 1 ) (59) 

oo 

= Y, e ~ TEn ^n(L2)MLi), (60) 

n=0 

where in the second line we used that the eigenstates form a complete orthonormal basis 
of the Hilbert space. Upon inserting ()53|) into ()60|) one obtains the unmarked finite time 
continuum propagator (cf. Appendix ID|) 



AV ; V L X L 2 S inh(v / 2AT) \ smh(V2AT) J 



(61) 



where Ii(x) denotes the modified Bessel function of the first kind. The corresponding 
propagator with a marking on the initial (final) spatial boundary can then be obtained by 
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multiplying (fHTj) with L\ (L 2 ). The finite time propagator (|HT|) of two-dimensional CDT 
was first obtained in jjj. 

A remarkable fact to mention at this point is that, up to some marking ambiguities, 
the finite time propagator (JBT]) agrees with the result of the propagator obtained from a 
continuum calculation in proper-time gauge of two-dimensional pure gravity pE2l . 

Since all continuum results presented so far were calculated in the Euclidean sector, 
the remaining task is to define an inverse Wick rotation to relate the Euclidean sector 
with the Lorentzian sector of the theory A natural proposal is to analytically continue 
the continuum proper time T i— > —iT. Under this prescription the effective quantum 
Hamiltonian corresponds to a unitary time evolution, 

U(T,T>) = e lHT ~ T \ with H(L, ±) = -L^ - 2 ± + 2AL. (62) 

Further, the Euclidean partition function (|58|) and the Euclidean finite time propagator 
(16 1|) can then be simply related to their corresponding Lorentzian expressions. For the 
case of the finite time propagator, the expression in the Lorentzian sector reads, 

G A (L U L 2 - T) = \\4^— 7= h ( 2 VZAL ^ (Lorentzian s.). (63) 

M ; VM2 ^sin( v / 2AT) V ism(V2AT)J 

Here one can clearly see that the propagator possesses oscillation modes in the time variable 
T. 

As a final remark, notice that the definition of the inverse Wick rotation is not a priori 
clear. Besides the definition given above, one might also consider the analytic continuation 
in the cosmological constant A i— > —iA, which would also correspond to a continuation 
a t i — ► —idt on the time-like cutoff. Unfortunately, it does not lead to a unitary theory, for 
which reason we prefer the analytic continuation on the proper time variable T i— > — iT as 
the physically sensible inverse Wick rotation. 



2.6 Physical observables 

Having calculated the effective quantum Hamiltonian, the partition function and the prop- 
agator, there remains the question: What is the physics behind those expressions? As one 
can see from (|17|). the classical Einstein equations in two dimensions are empty. Neverthe- 
less, we obtained a non-trivial expression for the effective quantum Hamiltonian. Therefore, 
whatever dynamics it describes will be purely quantum, without a physical non-trivial clas- 
sical limit. More precisely, we will see that it corresponds to quantum fluctuations of the 
spatial length of the "universe" . 

Interesting observables to calculate are the expectation values of the spatial length and 
higher moments. Since the only dimensionful constant in the model is the cosmological con- 
stant with dimension [A] = length -2 , all dimensionful quantities will appear in appropriate 
units of A. Using the expression of the wave function ([53)1 . one can calculate expectation 
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values of spatial length and all moments 

(L m ) n = (n\L m \n)= I dfi(L)i \ ! ,(L) L'" i '„(!). (01.) 



A general expression for (jMj) is presented in Appendix [Dl Let us just state the first two 
expressions here, namely, 

<L>„ = lt£, (L%=\<Z±£. (65) 



Generally, the moments scale as (L m ) n ~A 2 as shown in the appendix. From (|65|) one 
can also calculate the variance (AL) of the spatial length, 



(Ai>„ = \J(L 2 ) n - (L) 2 n = \ ?±±. (66) 



2 VA ' 

One can see that for the "universe" in a certain state | n), for example, the ground state 
I 0), it possesses fluctuations around the average value (L) with variance (AL) Q . This 
behavior is illustrated in Figure |H1 by a numerical Monte-Carlo simulation. What is also 
good to see in the figure is that the fluctuations are roughly of the same order as the spatial 
length of the "universe", i.e. (AL)~(L). 

Other useful quantities to calculate are certain critical exponents of the continuum 
theory, which, due to the gravitational setting, will all have a geometrical interpretation. 
Among those critical exponents, there is one of special interest, namely, the Hausdorff 
dimension d H (cf. |H]). In particular, dn can be calculated for each geometry [g^], but 
what we will consider here is the ensemble average over the entire class Geom(M) , that is, 
the leading order scaling behavior of the expectation value 

(V) T ~T dH . (67) 

To calculate this quantity we use the partition function (|58|). whose behavior for large T 
is given by 

Zr(A) ~ e ~ 2V1XT for T -> 00. (68) 
From this we can define a typical energy scale M(A) of the "universe" , 

M(A)= lim P° g2T(A) VAi (69) 



On the other hand, we can use the partition function Zr(A) to calculate the ensemble 
average of the volume, by taking derivatives with respect to the corresponding coupling, 
namely, the cosmological constant A, 
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Figure 6: A typical two-dimensional Lorentzian space-time. The compactified direction 
shows the spatial hypersurfaces of length (L) and the vertical axis labels proper time T. 
Technically, the picture was generated by a Monte-Carlo simulation, where a total volume 
of N = 18816 triangles and a total proper time of t = 168 steps was used. Further, initial 
and final boundary has been identified. 



Upon inserting (|68|) into (J7UJ), we can determine the scaling behavior of the volume at large 
scales T>1/M(A), 

<V> T ~T.A-4 for (71) 

This shows that at large scales T ^> 1/M(A), the typical universe has a volume which is 
proportional to T and therefore looks like a long tube of length T. Further, from (fTTj) one 
can read off that the spatial length scales as 

(V) 1 

(L) = -7^- ~ -=, 72 

as we have already found in ()65|). More interesting than the scaling of the volume at large 
scales is the typical scaling at small scales, namely, those of order T ~ 1/M(A). In the 
same procedure one obtains 

{v) ?~ r for T ~WAy (73) 



From definition (J67|) . one then sees that the Hausdorff dimension of 2d Lorentzian quantum 
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gravity or CDT is given by 

d H = 2 (CDT). (74) 

Naively this looks trivial, since we were considering two-dimensional CDT which in the 
construction uses two-dimensional building blocks. However, the Hausdorff dimension is 
truly a dynamical quantity and is a priori not the same as the dimension of the simplicial 
building blocks. An example of this is two-dimensional Euclidean quantum gravity 
defined through dynamical triangulations which has a Hausdorff dimension of dn = 4. 11 
This non-canonical value of the Hausdorff dimension is due to the highly degenerate, 
acausal geometries contributing to the ensemble average. At this point one can clearly 
see, as already mentioned in Section |2~3~1 that the continuum theories of two-dimensional 
quantum gravity with Euclidean and Lorentzian signature are distinct and cannot be 
related in a trivial manner (cf. footnote [7J) ! 

To conclude, we have seen that CDT gives a powerful method to define a nonperturba- 
tive path integral for quantum gravity. In the considered two-dimensional model with fixed 
cylindrical topology, it led to a well-defined continuum quantum theory, whose time evolu- 
tion was unitary. Further, we have seen that the restriction to the causal geometries in the 
path integral led to a physically interesting ground state of a quantum "universe" , which 
also showed that quantum gravity with Euclidean and Lorentzian signature are distinct 
theories. 

One can now ask the question whether or not causal constraints might be a helpful tool 
to implement topology changes and the sum over topologies in the path integral. Finding 
a detailed answer to this question will be the topic of the following Section EH 

Another aspect we have not considered so far is the question of universality of the 
theory: In general one would expect that the continuum theory should be independent of 
the details of the discretization procedure we use, e.g., whether we build up our space-time 
from triangles, squares or p-polygons should not affect any continuum physical observable. 
A discussion of this point will be the topic of Section HJ 



3 2D Lorentzian quantum gravity with topology 
change 

In this section we construct a combined path integral over geometries and topologies for 
two-dimensional Lorentzian quantum gravity. The section is structured as follows: In 
Section l3~T1 we introduce the reader to the general concepts of a path integral formulation for 
quantum gravity including a sum over topologies and discuss certain problems, especially 
in the context of Euclidean formulations. In Section I3~2l we demonstrate qualitatively how 
the Lorentzian structure can be used to exclude certain geometries from the path integral 
which lead to macroscopic causality violation; this makes the path integral well-behaved. In 

n One can rapidly check this by noting that the partition function in this case scales as Zt(A) ~ e - c \/at 
for T ^ oo |7|. 



3.1 Sum over topologies in quantum gravity 
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Section EP1 an explicit realization of this construction in the framework of two-dimensional 
CDT and the discrete solution thereof is presented. The difficulty in the construction of 
the continuum limit is to find a suitable double scaling limit for both the cosmological and 
Newton's constant, which we address in Section 13.41 The resulting continuum theory of 
quantum gravity describes a quantum "universe" with fluctuating geometry and topology, 
as we will discuss in Section 13.51 in terms of its physical observables. Further, we show in 
Section I3~H1 that the presence of infinitesimal wormholes in space-time leads to a decrease 
in the effective cosmological constant. 

3.1 Sum over topologies in quantum gravity 

In Section |21 we have given a concrete meaning to the nonperturbative gravitational path 
integral as an integral over all possible geometries Geom(M), 

Z= [ V[g^]e lS ^\ (75) 

JGeom(M) 

which by virtue of the approach of causal dynamical triangulations led to a well defined 
continuum theory of quantum gravity in two dimensions. 

Nevertheless, many other attempts of constructing a nonperturbative gravitational path 
integral start from the ansatz that one should not only consider all geometries of a manifold 
M of a certain fixed topology, but further include the sum over all space-time topologies 
in the gravitational path integral, formally written as 

Z = Y.\ fW^ W - (76) 

M JGeom(M) 

But how should we think about summing over topologies and further about the associated 
topology changes [T7] of space as a function of time? Clearly, theories which predict topol- 
ogy changes at macroscopic scales are unlikely to be consistent with observational data. 
A suggestion for how one should think about topology changes is therefore as topological 
excitations of space-time at very short scales such as the Planck scale. This leads to the 
notion of space-time foam )4*3) I44j . according to which space-time is a smooth manifold at 
macroscopic scales, but at small scales is dominated by roughly fluctuating geometries and 
topologies. Whether nature really admits this behavior and one should include the sum 
over topologies into the path integral is still unclear. 

The gravitational path integral including the sum over topologies has been the subject 
of intensive study in the context of four-dimensional Euclidean quantum gravity [JSJ ES] . 
However, in nonperturbative quantum gravity models, where one can analyze the problem 
explicitly, it turns out that topology changing configurations completely dominate the path 
integral 1)760 . since the number of contributing geometries grows super- exponentially with 
the volume N, i.e. at least ~iV!, which makes the path integral badly divergent. 

Moreover, this problem is still present in Euclidean quantum gravity in dimensions 
d<4. In the case of two dimensions, Euclidean quantum gravity without topology changes 
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is well understood analytically, as mentioned in Section El in the context of dynamical 
triangulations. Furthermore, the sum over topologies takes a very simple form, as a sum 
over a single parameter, 0^0, the genus of the space-time manifold M, i.e. the number of 
space-time handles. Including the sum over topologies, the Euclidean analogue of (|75|l has 
been studied intensively, also because it can be interpreted as a nonperturbative sum over 
worldsheets of a bosonic string with zero-dimensional target space (Uj. There it turns out 
that the divergence of the path integral still persists, moreover, the topological expansion 
of ()7fij) is not even Borel-summable |18j . Further attempts to solve this problem in the 
Euclidean context have so far been unsuccessful. 

To draw a lesson from Section |2l we have seen that in pure Lorentzian quantum gravity 
or CDT (we will from now on refer to the CDT model with fixed topology as the "pure" 
model) the restriction to causal geometries in the gravitational path integral (J75|) leads to 
a physically interesting ground state of a quantum "universe" which is distinct from the 
Euclidean case. Inspired by this, one might now ask the question whether implementing 
the (almost everywhere) causal structure into the gravitational path integral including 
the sum over topologies might restrict the state space in such a way that the number of 
geometries does not grow super-exponentially anymore. That this is indeed the case, and 
that one can unambiguously define ()76)1 in the setting of CDT has been shown in |2*T] 12*2*] 12 
and will be discussed in the following sections. 

3.2 Implementing topology changes 

Having adopted the notion of space-time foam and the sum over topologies in the path 
integral, we have to face the question, what kind of topology changing geometries should 
contribute to the sum in (|7fi|). 

From observational data we tend to exclude all topology changing geometries which 
lead directly to a macroscopic acausal behavior. In our implementation we therefore con- 
sider only topology changes which are of an infinitesimal duration, so-called infinitesimal 
wormholes. In terms of the framework of CDT this means that wormholes do only exist 
within one time step At = 1 of discrete proper time. Further, the total number of worm- 
holes per such discrete time step can be arbitrary in the continuum limit. Using this notion 
of infinitesimal wormholes, we consider triangulations whose spatial slices have topology 
S* 1 for integer values of t, whereas in the interval }t, t+l[ it splits into flt+1 ^-components, 
where g t is the number of wormholes or equivalently the genus of this space-time strip. 

The setting with infinitesimal wormholes in a path integral formulation of quantum 
gravity with a sum over topologies has been studied before in the context of Euclidean 
quantum gravity in four dimensions fI5\ |2*E] . Nevertheless, as we will show in the fol- 
lowing, even generic space-times adopting this mildest form of topology change can cause 
a super-exponential growth in the number of configurations in the path integral or lead 
to a macroscopic acausal structure. The essential difference between Euclidean quantum 
gravity and CDT is the (almost everywhere) Lorentzian structure, which enables us to 



12 Also see [231 IS] for additional information. 
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Figure 7: (al) A spatial slice of topology S 1 at time t splits into three components giving 
rise to two saddle points p t . (a2) A merging of two components at time t + l, where the 
saddle points p t and p t+ \ are time-like related. This gives rise to a "untwisted" wormhole, 
as shown in (a3). Parallel light rays passing this wormhole are unaffected unless they are 
not scattered directly by the wormhole, as emphasized in (a4). For the case of a twisted 
regluing (bl), where the saddle points p t and p t+ i are space-like separated, parallel light 
beams passing this wormhole split into two parts with a relative separation (b2). 



classify how badly causality is violated and to exclude certain "badly" causality violating 
geometries from the path integral. To be more precise on this point, clearly the metric 
of the topology changing geometries becomes degenerate at the saddle points p, where 
space-time splits or merges. Nevertheless, one is still able to assign a light-cone to these 
so-called Morse points by analysis of the neighboring points as discussed in jlHl HH EH EI] ■ 
Let us now give some qualitative arguments for excluding certain classes of geometries 
with "bad" topology changes, and then present the concrete realization of how to implement 
the remaining "not-so-bad" topology changes in the framework of CDT in the next section. 

The creation of wormholes is illustrated in Figure At discrete proper time t the 
spatial slice of topology S 1 splits into flt + 1 S^-components giving rise to Q t saddle points 
(al). After a time lapse At — 1, the S^-components merge again to a single S 1 , which leads 
to Qt wormholes in the space-time strip [t, t+l]. An important point in this construction is 
the question of how one is going to reglue the different ^-components at time t+l. It is not 
difficult to see that arbitrary regluings at time t + l give again rise to a super-exponential 
growth of the number of configurations. However, most of these regluings are very ill- 
behaved with respect to their "causal structure" , in the sense that part of light beams can 
get a global rearrangement when passing through the wormholes. Fortunately, the imposed 
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(almost everywhere) Lorentzian structure enables us to restrict ourselves to a subclass of 
topology changes which does not admit this behavior, namely, regluings without a relative 
rearrangement or twist of the corresponding S^-components. These topology changes are 
characterized by the fact that the upper saddle point p t +i is time- or light-like related to 
the lower saddle point p t , as indicated in Figure [3 (a2)-(a3), where for simplicity only one 
merging of two S^-components is shown. 

To illustrate the qualitative difference between those untwisted and twisted topology 
changes, consider light beams propagating through the resulting space-times, as shown in 
Figure [?l (a4) and (b2). In both cases light beams which hit the wormhole scatter non- 
trivially. However, in the twisted case another effect takes place, namely, due to the relative 
spatial shift Al between the saddle points p t and Pt+i, parallel light beams passing this 
wormhole will split into two parts with a relative separation AL While the effect due to 
the non-trivial light scattering goes to zero in the continuum limit, the relative separation 
in the twisted case persists, even after the wormhole has disappeared. We therefore discard 
these twisted mergings as "bad" topology changes and exclude them from the path integral. 
Thus, we only consider geometries with untwisted wormholes in the sum over topologies. 
Note that even though the effect of light scattering at such an untwisted wormhole goes to 
zero in the continuum limit, this does not mean that the contribution of several wormholes 
will not lead to a measurable effect; in fact, this is precisely what we will observe in the 
resulting theory. 

3.3 Discrete solution: The transfer matrix 

In this section we discuss how wormholes of the type explained above can be introduced 
in the setting of CDT. Further, we present the combinatorial solution of the one-step 
propagator including the sum over topologies in the discrete setting. 

Consider a typical space-time strip of topology [0, ljxS" 1 , as used in the discrete setting 
of the pure CDT model, which is illustrated in Figure |H1 (i)- An "untwisted" wormhole of 
infinitesimal duration At = l, as associated to a "not-so-bad" geometry, can be created by 
identifying two time-like edges of this space-time strip (ii) and cutting open the geometry 
perpendicular to this line (iii). Clearly, the resulting saddle points p t and p t +i are time- like 
related. The curvature singularities at those saddle points will be of the standard conical 
type after we have performed the Wick rotation and we can assign the Boltzmann weight 
accordingly. 13 A space-time strip of arbitrary genus q can then be generated by repeating 
this procedure g-times, where the q arrows identifying time-like edges are not allowed to 
intersect, as illustrated in Figure |H1 This type of identification ensures that we obtain just 
an exponential growth in the volume, which leads to a well-defined continuum theory, in 
contrast to the super-exponential growth in the Euclidean model. 

As in the pure CDT model, all dynamical information is encoded in the transfer matrix 
and it therefore suffices to investigate the combinatorics of a single space-time strip. Hence, 
following the prescriptions of Section 12.41 and including the topological term (j!7|) in the 
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Figure 8: Construction of a wormhole: starting from a space-time strip of topology 
[0, 1] x 5* 1 as in the pure CDT model (i), one identifies two time-like edged (ii) and then 
cuts open the geometry perpendicular to this line (iii). The two resulting saddle points at 
time t and t+1 are labeled with p t and pt+i- 



Regge action, we obtain the Wick rotated one-step propagator 

TeW(T): 
'in ^out 

where the sum is taken over all possible triangulations T of height At = 1 with fixed 
initial boundary l in and final boundary l out , but arbitrary genus ^ &{T) ^ [N/2] (here 
N = li n + lout is the number of triangles in T, which coincides with the number of time- 
like edges). Further, k is the bare inverse Newton's constant and A the bare cosmological 
constant. As in the pure model we have omitted an overall constant phase factor. 

Let us now perform the sum in (J77j) for fixed genus. Consider first the case without 
wormholes. To simplify the combinatorial expressions, we will use the combinatorial factor 
belonging to a space-time strip with open boundary conditions, instead of circular boundary 
conditions. As shown in Appendix [HIP the resulting continuum theories in the pure CDT 
model are very similar for both boundary conditions, which enables us to recover the result 
for circular boundary conditions at a later stage. Hence, the combinatorial factor assigned 
to the sum over all possible triangulations without wormholes reads (Appendix [BJ 

G(kn, Lut) = ( m j ° J ■ (78) 

\ 'in / 

For a given strip with fixed number of triangles N = li n +l ou t, wormholes are created by the 
procedure explained above. Thereby, to construct a space-time strip with 0^$j(T) ^ [N/2] 
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appear as dots on a circle in the cutting plane. Pairwise identification of the 2g time-like 
edges (dots) is illustrated by the non-intersecting arches, where maximally one arch can 
connect to each dot. 



wormholes, one first has to choose 2g out of iV time-like edges; the combinatorial factor 
assigned to this is simply given by 

For a given set of 2q time-like edges we have to count the number of possibilities to pairwise 
identify these edges, where the arches belonging to the identifications are not allowed to 
intersect, as illustrated in Figure |H1 This is a well-known combinatorial problem whose 
solution is given by the Catalan numbers 



(20)! 1 f2g 



fl!(fl + l)! 0+1 W 



SO) 



Thus, the complete formula for the one-step propagator (|77|). after the fixed genus part of 
the sum has been performed, reads 

From this one can obtain the propagator G\ tK (li n , l out ; t) for arbitrary time t by iterating 
(18 1|) t times according to the composition law 

G\ lK (li n , lout: ti + t-z) = Gx >K {kn, I, t\)G\ >K (l, lout, ^2); (82) 

I 

Gx^lmjloutjt + 1) = G\ >K (lj n , I, l)G\ tK (l, Iput, t). (83) 

I 

To give a complete solution of the discrete problem we still have to explicitly perform the 
sum over all genus in (|81|). We will obtain this result in Laplace transformed "momentum" 
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space, by introducing generating functions for the one-step propagator 

oo 

G{x,y;g,h;l) = G xA l in,loutA)x kn y lou \ (84) 

I in i^-out —0 

where we have defined g = e~ A , as in the pure model, and h = e~ K . Upon inserting (J%T|) 
into (J%%|) an d evaluating the summations over the Z's, one obtains the generating function 
of the one-step propagator 

1 2 

G(x, y, g, h, 1) = (85) 

1-5 + y) 1 + V l - 4m 2 

with 

u = - j . (86) 

Note that in order to arrive at this expression, we have explicitly performed the sum over 
all (not-so-bad) topologies! The fact that this infinite sum converges for appropriate values 
of the bare couplings has to do with the restriction to the "not-so-bad" topology changes 
according to the imposed causality constraints, as motivated in the previous section. 
In (185)1 one recognizes the generating function Cat(ii 2 ) for the Catalan numbers, 



„ / 2 x 2 1 - VI - 4w 2 . , 

Cat (it ) = - = V — 2 . 87 

1 + V 1 - 4m 2 2u z 

For h = one has Cat(w 2 ) = 1 and expression (|85j) reduces to the generating function of the 
one-step propagator with open spatial boundary conditions of the pure CDT model (cf. 
Appendix lB| . 

G(x, y;g,h = 0; 1) = . (88) 

1 - gx - gy 

3.4 Continuum and double scaling limit 

Taking the continuum limit in the case of pure CDT is fairly straightforward, as we have 
seen in Section 12.51 The joint region of convergence of the one-step propagator (|88|) is 
given by 

\x\ < 1, \y\ < 1, \g\ < -. (89) 

The continuum limit was then obtained by simultaneously fine-tuning to the critical values 
with use of the following canonical scaling relations 

g = I(l-a 2 A) + 0(a 3 ), (90) 

x = 1- aX + 0(a 2 ), (91) 
y = 1-aY + 0(a 2 ) } (92) 
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where X, Y and A denoted the renormalized couplings. The difficulty which arises when 
taking the continuum limit in the case with topology changes is to find a suitable double 
scaling limit for both the gravitational and cosmological coupling, which leads to a phys- 
ically sensible continuum theory. By this we mean that the one-step propagator should 
yield the Dirac delta-function to lowest order in a, and that the Hamiltonian should be 
bounded below and not depend on higher-order terms of the scaling relations for all cou- 
plings, in a way that would introduce a non-trivial dependence on new couplings in the 
continuum theory. Difficulties in finding such a double scaling limit arise due to the fact 
that Newton's constant is dimensionless in two dimensions, whence there is no preferred 
canonical scaling for h. One can make the multiplicative ansatz 14 

h = —=h ren (ad) p , (93) 
V2 

where h ren depends on the renormalized Newton's constant Gn according to 

In order to compensate the powers of the cut-off a in (}9~3*]) . d must have dimensions of 
inverse length. The most natural ansatz in terms of the dimensionful quantities available 
is 

d= (VX a (X + Yf- a ) } (95) 

where the constants a and (3 have to be chosen such that one obtains a physical sensible 
continuum theory according to the above considerations. 

To calculate the effective quantum Hamiltonian one can follow a similar procedure as 
the one used in the case of the pure CDT model in Section 12*31 Thereby, we start off with 
the one-step time evolution of the discrete wave function, recall 



*l>(x, t + l)=<f —^G(x, z- 1 ; g, h; 1) ^(z, t). (96) 
J lixiz 

Upon inserting the scaling relations ()90 |1 -(|92 |1 and t = — into this equation and using 

ip{X, T + a) = e- aA{x ' 9x) ij(X, T), (97) 
one can expand both sides to first order in a, yielding 

(l - *H + £,(„>)) , K X) = £ |§ { + a Cat M } ,(Z), (98) 

with ip{X) = ip(x = l — aX). For convenience, we treated separately the first factor in the 
one-step propagator (|55)i. which is nothing but the expansion of the one-step propagator 



14 Here the factor is chosen to give a proper parametrization of the number of holes in terms of 
Newton's constant (see Section ET5)l . 
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without topology changes, and the second factor, the Catalan generating function (|S7jl. 
which contains all information on the new couplings. Note that the first term on the right- 
hand side of ()98|1. -^x^ * s ^ ne Laplace-transformed delta-function. This gives us sensible 
information of the values of a and f3 in the expansion of the Catalan generating function. 
Inserting the scaling relations (j9*Uj) - (}9*3"j) into Cat(w 2 ) and expanding, one obtains 

Cat( M 2 ) = 1 + 2 ^ k ^ en 2 a w ~ 2 + 0{a 2f3 - 1 ). (99) 



Thus, in order to preserve the delta-function to lowest order in (J98|) and to have a non- 
vanishing contribution to the Hamiltonian one is naturally led to (3 = 3/2. For suitable 
choices of a it is also possible to obtain the delta-function by setting (3 = 1, but the 
resulting Hamiltonians turn out to be unphysical or at least do not have an interpretation 
as gravitational models with wormholes, as we will discuss in Appendix IE1 15 

Considering /3 = 3/2 and inserting ((55]) into (|5S|). the right hand side of this equation 
becomes 

ix 2ni\z-x \(x-zf (x-zf)y y ' { ' 

Observe that to have a non-trivial dependence on the new couplings in the continuum 
theory, one should not allow for values a^O. Hence, performing the integration for a>0 
and discarding possible fractional poles in (jlOOj) . the effective quantum Hamiltonian in 
"momentum" space reads 

^•W^ x ^ + x - 2A W + 2A ^Tk^- "4 1' 1 -" (101) 

For all possible values of a, these Hamiltonians do not depend on higher order terms in the 
scaling relations (|50 ]l -(p3 ]l . One might check this, as explained in Section 1231 by explicitly 
introducing higher order terms in the scaling relations and observing that the resulting 
Hamiltonians are independent of them. 

One can now perform the inverse Laplace transformation on the wave function, 
ip(L) = 7^.e x L ip(X), to obtain the effective quantum Hamiltonian in length space, 

8 8 2 8 Itsr^h 2 

H ^M^- L m?-TL + 2AL -T^ L ''' —M,2,..., (102, 

where we defined m = 3a— 1. Since H is unbounded below for m^2, we are left with m = 
and m = l as possible choices for the scaling. However, setting m = merely has the effect 
of adding a constant term to the Hamiltonian, leading to a trivial phase factor for the wave 
function. Hence, we conclude that the only continuum theory with non-trivial dependence 
on the new couplings corresponds to m = l, i.e. to the following double scaling limit 

fc* = \h 2 ren A (X + Y) a 3 , (103) 



15 One might also consider scalings of the form h — * cih ren (ad)+C2h ren (ad) 3 ^ 2 , but they can be discarded 
by arguments similar to those of Appendix El 
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with the Hamiltonian given by 

H(L, A) = - L A_ _ JL + (i _ 2 A L (1Q4) 

Note that for all values Gn > of the renormalized Newton's constant ()93|) the Hamiltonian 
is bounded from below and therefore well defined. Further, the Hamilton operator is 
selfadjoint on the Hilbert space TC = £ 2 (M + , dL). Note that (jl(J4|) corresponds to the 
propagator with open boundary conditions. It is not difficult to recover the Hamilton 
operator for the case of circular boundary conditions 

dL ] = ~ L dV ~ 2 dL + ^ " k2 ™) 2AL - ( topology S ' X t°' ^ ( 105 ) 

which is selfadjoint on the Hilbert space 7i = £ 2 (IR + , LdL). The spectrum of this Hamil- 
tonian is discrete with equidistant eigenvalues, 

E n = 2v/2A(l-/i r 2 en )(n + l) 71 = 0,1,2,.. (106) 

The corresponding eigenfunctions are given by 



ML) = Xe'^d-^Ji iFi( _ n . 2 . 2v /2A(l-^ en )L), = LdL, (107) 

where i-Fi(— n, fi + 1; z) are the confluent hypergeometric functions, defined in (|55j). As in 
the case of pure CDT, the eigenfunctions ip n {L) form an orthonormal basis of the Hilbert 
space Ti, where the normalization factor is given by 



An = 2v/2A(l-/4J(n + l). (108) 

Having obtained the spectrum of the Hamiltonian one can easily calculate the (Eu- 
clidean) partition function 



ZHA.Gd-E'-"- = ^=J= F - , ^ = e~>°'. (109) 
Further the time propagator or loop-loop correlator can be obtained as (cf. Appendix |D|) 

oo 

Ga,g n (Li,L 2 ;T) = J2 e ~ TE "^n( L 2)MLi) (110) 



n=0 

u e -a,(L 1+ L 2 )coth(a;r) / 2 Uy/LJ^\ 

y/Tj7 2 sinh(wT) 1 V sinh(a;T) J ' ^ n ) 

where we have used the short hand notation uo = a/2A(1 — /i 2 en ). As expected, for /i ren — > 
the results reduce to those of the pure two-dimensional CDT model, as obtained in Section 
1231 
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3.5 Physical observables 

In the following we want to analyze the new model of CDT with topology changes in 
terms of its physical observables. As in the pure model, there are certain geometric 
observables such as the expectation values of the spatial length and higher moments, 
(L m ) n = (n \ L m \ n), which can be calculated using the expressions for the wave functions 
(jl()7j) . as obtained in the previous section. The first results read 

~ a/2A(1 - h 2 ren ) ' (L) "-4A(l-^ en )' (112) 

where a general expression for higher moments is presented in Appendix El Generally, the 
moments scale as (L m ) n = A~^" as in the pure CDT model. From (jll2|) one can then also 
calculate the variance (AL) n of the spatial length, 

(AL) »=^^^^ra (113) 

One can still interpret these observables in terms of a fluctuating universe, as done in the 
pure CDT model, where the quantum fluctuations are now governed by the "effective" 
cosmological constant A cff = A(l — h^ en ). A detailed discussion on this point in terms of 
some new "topological" observables will be given in the next section. 

Another geometrical observable is the Hausdorff dimension djj, which can be calculated 
from the ensemble average over all geometries with arbitrary genus. Since the partition 
function scales as 

Z T (A, G N ) ~ eV 2A ( 1 - ft ™) T as T — > oo, (114) 

a similar calculation, like the one performed in Section I2.6| reveals that the Hausdorff 
dimension is given by 

d H = 2. (115) 

This is the same result as obtained in the case of pure two-dimensional CDT. 

In addition to these well-known geometric observables, the system possesses a new type 
of "topological" observable which involves the number of wormholes N s , as already antici- 
pated in Section I3.2[ where we mentioned that the presence of wormholes in the quantum 
geometry and their density can be determined from light scattering. An interesting quan- 
tity to calculate is the average number of wormholes in a piece of spacetime of duration 
T, with initial and final spatial boundaries identified. Because of the simple dependence 
of the action on the genus this is easily computed by taking the derivative of the partition 
function Zt(A,Gn) with respect to the corresponding coupling, namely, 

/ a j v 1 Ken OZ T (A, G N ) 

™ Z T (A,G N ) 2 dh ren ■ 
Upon inserting (jl(J9J) this yields 

coth( v /2A(l-/ i 2 en )T)-l 
(N g ) = T h 2 ren A V - =J- . (117) 
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Figure 10: The density of wormholes n in units of A as a function of Newton's constant 
Gn- 



In an analogous manner we can also calculate the average spacetime volume 



(V)- 1 dZ T (A,G N ) 

{V) -~Z T (A,G N ) M ' (118) 



leading to 



( V) = T y 1 2A reJ (coth (^AMJT) - l) . (119) 
Dividing ()117j) by ()119j) we find that the spacetime density of wormholes n is finite, 

"=w = i^: A - (120) 

The density of holes in terms of the renormalized Newton's constant is given by 

1 

e G N — 1 



n=— A. (121) 



The behaviour of n in terms of the renormalized Newton's constant is shown in Figure 
ITUl The density of holes vanishes as Gn — »0 and the model reduces to the case without 
topology change. 



3.6 Taming the cosmological constant 

We have already seen in the last section that fluctuations in the geometry gave rise to 
the notion of the "effective" cosmological constant A e s = A(l — h% en ). In the following we 
want to interpret this effect in terms of the physical quantities, namely, the cosmological 
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Figure 11: The coefficient of the effective potential, f(rj) = 1/(1 + 7/), as function of the 
density of holes in units of A, rj — j. 



scale A and the density of wormholes in units of A, i.e. V — T- These can be seen as the 
two scales of the model, in contrast to the pure CDT model which just has a single scale. 
As in the latter, the cosmological constant, as the only dimensionful quantity, defines the 
global length scale of the two-dimensional "universe" through (L) ~ 4=. The new scale 
in the model with topology change is the relative scale rj between the cosmological and 
topological fluctuations, which is parametrized by Newton's constant Gn- Both together 
govern the effective fluctuations in the quantum geometry. 

One can nicely observe this when rewriting and interpreting the Hamiltonian (jl()4j) in 
terms of the new physical quantities, A and 77, resulting in 

d d 2 d 1 

h ^m> = - l m?-Tl + —, 2KL - < 122 > 

One sees explicitly that the topology fluctuations affect the dynamics since the effective 
potential depends on 7/, as illustrated by Figure [TTJ From the coefficient of the effective 
potential, f(rj) = 1/(1—7/), one also observes that the presence of wormholes in our model 
leads to a decrease of the "effective" cosmological constant A e ff = /(?/) A. 

This connects nicely to former attempts to devise a mechanism, the so-called Coleman's 
mechanism, to explain the smallness of the cosmological constant in the Euclidean path 
integral formulation of /cmr-dimensional quantum gravity in the continuum with the pres- 
ence of infinitesimal wormholes [2H1I2S1- The wormholes considered there resemble those 
of our toy model in that both are non-local identifications of the spacetime geometry of 
infinitesimal size. The counting of our wormholes is of course different since we are working 
in a genuinely Lorentzian setup where certain causality conditions have to be fulfilled. This 
enables us to do the sum over topologies explicitly. 

Further, in Coleman's mechanism for driving the cosmological constant A to zero, 
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an additional sum over different baby universes is performed in the path integral, 
which leads to a distribution of the cosmological constant that is peaked near zero. 
We do not consider such an additional sum over baby universes, but instead have 
an explicit expression for the effective potential which shows that an increase in the 
number of wormholes is accompanied by a decrease of the "effective" cosmological constant. 

To conclude, we have seen that two-dimensional Lorentzian quantum gravity including a 
sum over topologies leads to a well-defined unitary continuum quantum theory, different to 
the one of the pure CDT model. The presence of causality constraints imposed on the path- 
integral histories - physically motivated in Section enabled us to derive a new class of 
continuum theories by taking an unambiguously defined double-scaling limit of a statistical 
model of simplicially regularized space-times. These causality constraints were crucial to 
resolve the problem of super-exponential growth in the number of configurations, as present 
in the model of two-dimensional Euclidean quantum gravity with topology changes. The 
resulting model of two-dimensional CDT with topology changes has, besides the well- 
known geometrical observables, new "topological" observables, such as the finite space- 
time density of wormholes. Further, we observed that the presence of wormholes in our 
model leads to a decrease of the "effective" cosmological constant, which nicely connected 
to former attempts of driving the cosmological constant A to zero in the formulation of 
/ottr-dimensional quantum gravity in the continuum through the presence of infinitesimal 
wormholes, as described by Coleman's mechanism. 



4 Universality: Its physical and mathematical impli- 
cations 

In this section we discuss generalizations and implications of the "pure" two-dimensional 
CDT model, as described in Sectional In doing this, we use the techniques of the formalism 
developed in Section |2l The section is structured as follows: In Section l4~Tl we introduce a 
higher curvature term in the pure model and show that the resulting theory is equivalent to 
a two-dimensional Lorentzian quantum gravity model whose discrete space-times consist of 
squares and triangles. Further we show that it belongs to the universality class of the pure 
two-dimensional CDT model. In Section 14.21 a model is introduced which only allows for 
minimal curvature weights per space-time strip, where we see that the resulting model also 
shares the same universality class with the pure two-dimensional CDT model. This makes 
it clear that two-dimensional CDT define a rather broad universality class and is thus a 
genuine continuum theory independent of the details of the regularization. In Section |4~3*1 
we discuss bijections between CDTs, heaps of dimers and Dyck paths. Such connections 
can be a useful mathematical tool to reduce the two-dimensional counting problem to an 
effective one-dimensional problem. 



4.1 2D Lorentzian quantum gravity with a higher curvature term 
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(a) (6) 

Figure 12: (a) A typical triangulation with staircase boundary conditions. The dashed 
lines show the relation to the corresponding dual triangulation which is drawn in (b) . 



4.1 2D Lorentzian quantum gravity with a higher curvature term 

Following [53] we introduce a higher curvature term in the pure model of 2D Lorentzian 
quantum gravity. We will see that the resulting theory preserves universality of the pure 
model. Further, the higher curvature term can be interpreted as being equivalent to in- 
serting squares in the triangulation of the pure model. 

In the pure two-dimensional CDT model we attached a "cosmological" weight g = 
e to every triangle face resulting to the factor 

gN(T) in 

the Boltzmann weight of each 
triangulation T, where N(T) = li n +l ut is the number of triangles in the triangulation. Since 
the usual curvature term in the Einstein-Hilbert action is trivial in 2d for fixed topology, 
there was no curvature weight in the propagator. 

We now want to generalize the pure model by explicitly introducing a higher curvature 
weight into the propagator which suppresses or enhances local curvature. To still be able 
to use the powerful techniques of the transfer matrix formalism, we require that this higher 
curvature weight is defined locally within one strip of height At — 1. A convenient way to 
define such a higher curvature weight is by assigning a weight G to each pair of adjacent 
triangles pointing both up or pointing both down [SHI El]- The resulting effect of this higher 
curvature term in the eigenvalues of the transfer matrix can most easily be calculated 
using the dual graphical representation of the triangulation. For convenience we consider 
a triangulation with "staircase" boundary conditions 16 , as illustrated in Figure IT2*1 (a). In 
the dual representation (b) a space-time strip of height At = 1 translates into a sequence of 
half-edges attached to the dual constant time line, where the half-edges above (below) this 
line correspond to the up-pointing (down-pointing) triangles in the triangulation (a). In 
terms of the dual triangulation we assign a cosmological weight g to every vertex (which is 



16 Thc resulting theory for "staircase" boundary conditions is closely related to the one with periodic 
boundary conditions as shown in Appendix iBl 
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in fact the dual of the triangle face) and in addition each pair of neighboring up-pointing 
half-edges or down-pointing half-edges receives a factor 0, 



1 111 . (123) 

gg i i e e 

Using these weights we can calculate the eigenvalues of the transfer matrix 
G g ,e(hnJouut — 1) = {lout I T | l in ) by summing over possible configurations of a dual 
time line with l in half-edges below and l out half-edges above. In spirit of the higher cur- 
vature weight this is done by first summing over the number k ^ of blocks of n r ^ 1 
neighboring lower half-edges and m r ^ 1 neighboring upper half edges, with r = l,...k, 
and summing over all partitions of lower half-edges, U n = ^2 r n r , and upper half-edges , 
hut = J2r m ri yielding 

G g , e (hn, hut; t=l) = (g Q)^ £ ^_ (^~^ . (124) 

Introducing the generating function on the propagator G gt Q(l in , l out ; t = l) one obtains 
G(x,y;g,Q;l) = ^ G g,e(hn,lovt',t=l) 



g 2 xy 



1 — g Q(x + y) — g 2 (l — Q 2 )xy 



(125) 



As expected, for = 1 one recovers the one-step propagator of the pure model (jB.8|) . 
whereas for = curvature is fully suppressed and one gets the one-step propagator for a 
flat triangulation. 

It is very interesting to see that the discrete one step-propagator (|125j) can be reinter- 
preted as the one-step propagator of a CDT model whose discrete space-times consists of 
triangles and squares. 

Expanding the denominator of ()125j) . one only gets sequences of the terms gOx, gQy 
and g 2 (l — O 2 ). We can now assign a factor g<dx (gQy) to a lower (upper) half-edge which 
both have a weight gO. The factor g 2 (l — O 2 ) then belongs to a pair of upper and lower 
half-edge which can be merged to a crossing. The resulting new weights read, 



n • ( 126 ) 

g e g 2 (i-e 2 ) 



The corresponding space-time is shown in Figure EI] (a). The dual of this space-time is 
made up of triangles and squares, as illustrated in (b). In the following we want to look at 



4.1 2D Lorentzian quantum gravity with a higher curvature term 
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(a) (6) 

Figure 13: Reinterpretation of the resulting model with higher curvature terms. We assign 
a weight gQ to isolated half-edges, _L and T, and a weight g 2 {l — 2 ) to crossings +. A 
typical triangulation corresponding to these weights is drawn in (a). The dual triangulation 
(b) consists of up- and down-pointing triangles and squares. 



the continuum theory of such a two-dimensional Lorentzian quantum gravity model whose 
space-times are made up of triangles and squares in the discrete setting. 

Performing the continuum limit of ()125j) is straightforward. The joint region of conver- 
gence of (|125j) is given by 

|x| < 1, \y\ < 1, \g\ < -j-q-Q- (127) 

One can fine-tune to the critical point (x c ,y c ,g c ) = (1,1, j+q) by use of the following 
scaling relations 

g = -±^(l-a 2 A) + CV), (128) 

x = 1- aX + 0(a 2 ), (129) 
y = 1- aY + 0(a 2 ). (130) 

Upon inserting (J128j) - (J13()j) into (|125|) and following the formalism developed in Section El 
one obtains the effective quantum Hamiltonian of the resulting continuum theory, 

*( L 4» = - ei |j- e il + 2Ai - (131) 

which is self-adjoint on the Hilbert space TC = C 2 (R + , dL). Further, using (|D.28jl . one can 
read off the finite time propagator 

G A ,e A,L 2 ,T = J~ = J -== ^— ■ (132) 

' Ve sinh(VeAT) VVesinh {VqKt)J y ' 
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One sees that the resulting higher curvature model is equivalent to the pure CDT model 
(with staircase boundary conditions) up to a global rescaling of the coordinates, A\— > ^, 
Li^QLi and L 2 ^QL 2 . Hence, the new discretization of the two-dimensional Lorentzian 
quantum gravity model whose space-times are made up out of triangles and squares belongs 
to the same universality class as the pure CDT model which is discretized with triangles 
only. 

In jSHl h is shown that one can generalize this treatment to discrete space-times which 
incorporate p-polygons, where the resulting continuum theory still belongs to the same 
universality class as the pure two-dimensional CDT model. 

4.2 2D Lorentzian quantum gravity with minimal curvature 



In this section we want to discuss another generalization of the pure two-dimensional 
CDT model, where one only allows for minimal (positive or negative) curvature weights 
per space-time strip of height At = 1. By this we mean that the number of neighboring 
triangles pointing in the same direction is not allowed to exceed two. 

In terms of the discrete setting the construction is very simple: to calculate the eigen- 
values of the transfer matrix, G\(li n , l out ; t— 1) = (l out | T | l in ), one sums over all possible 
triangulations of a space-time strip of height At = 1 with l in up-pointing triangles and l out 
down-pointing triangles, where in addition one restricts the maximal number of neighbor- 
ing triangles pointing in the same direction to N max = 2. In the continuum limit this means 
that we only consider an infinitesimal amount of curvature per infinitesimal space-time 
strip. However, we expect that in a small, but finite region of space-time we can still 
generate arbitrarily large values of curvature. More precisely, we expect that the result- 
ing continuum theory will be equivalent to the one obtained in Section 14.11 with a fixed 
curvature suppressing factor Oe (0, 1). 

To solve the discrete problem we introduce the generating function for the one-step 
propagator G\(li n , l out ; t — 1). Hereby one can take advantage of graphical methods, where 
one assigns a factor of gx to every up pointing triangle "A" and a factor of gy to every down 
pointing triangle "y" • The one-step propagator for the case of open boundary conditions 
then reads 



weights 





k=0 v 

1 



1 - g 2 xy(l +gx)(l +gy)' 



(133) 



The joint region of conver; 



■g' 



;ence of (|133|) is given by 




(134) 



4.3 Causal dynamical triangulations as heaps of dimers and Dyck paths 



43 



To perform the continuum limit we fine-tune to the critical point 
(x c , y c , g c ) = (1, 1, \ {\f§ — 1)) by use of the following scaling relations 

9 = i(v / 5-l)(l-a 2 A) + 0(a 3 ), (135) 

x = 1-aX + 0(a 2 ), (136) 
y = 1 -aY + 0(a 2 ). (137) 

Inserting the scaling relations (jl35|) - (J137J) into the one-step propagator ()133|) and following 
the standard procedure one obtains the effective quantum Hamiltonian of the resulting 
continuum theory, 

^m ) = - eL &- e m +2AL ' withe =Srr a17 ' (138) 

which is self-adjoint on the Hilbert space TC = C 2 (R + , dL). One sees that, as in Section 
14. 1| the resulting model is equivalent to the pure CDT model (with staircase boundary 
conditions) up to a global rescaling of the coordinates, Ah| and L i— > QL and hence 
belongs to the same universality class as the pure CDT model. 

The above examples demonstrate that two-dimensional CDT define an entire univer- 
sality class, implying that it is independent of the subtleties of the regularization and is a 
genuine continuum theory of Lorentzian quantum gravity, without evidence of a discretiza- 
tion of space-time in a physical sense. 



4.3 Causal dynamical triangulations as heaps of dimers and Dyck 
paths 

In this section we want to discuss bijections between CDTs, heaps of dimers and Dyck paths 
which were first established in |)3] in a slightly modified form. Such connections can be 
useful to reduce the two-dimensional counting problem to a one-dimensional combinatorial 
problem. 

In the "pure" two-dimensional CDT model, as discussed in Section we considered 
space-time triangulations to have the topology [0, lJxS* 1 of a cylinder. Further, for combi- 
natorial reasons, we introduced a marking on the initial spatial boundary. One can follow 
the right-most time-like edge connected to the marking to trace a unique path through 
the entire triangulation which can be described by a sequence of up-pointing triangles. 
Cutting open the geometry along this path one can establish a one-to-one correspondence 
between cylindrical triangulations with a marking on the initial spatial boundary and tri- 
angulations of topology [0, 1] x [0, 1] whose left-most triangles are all pointing upwards, 
but without restrictions on the right-most triangles (cf. Appendix [BJ). A triangulation 
with such boundary conditions and the corresponding dual triangulation are illustrated in 
Figure HU 

Starting from a dual triangulation with this type of boundary conditions (Figure El 
(b)) we can remove the space-like links and view the remaining time-like links as dimers. 
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Figure 14: (a) A triangulation of topology [0, 1] x [0, 1] with all left-most triangles pointing 
upwards and no restrictions on the right-most triangles. Triangulations with these bound- 
ary conditions are in one-to-one correspondence with triangulations of topology [0, 1] x S 1 
with a marking on the initial spatial boundary, (b) The corresponding dual triangulation. 



Putting a base line to the left and pushing all dimers to the left until they touch the 
base line or another dimer we get a so-called heap of dimers, as illustrated in Figure 
ITol (b). The dual of this heap of dimers is a tiling of diamonds fFigure ITB1 (a)), which 
one can obtain by extending the original triangulation to a vertex below and above and 
then removing all space-like links. The combinatorial structure of heaps of dimers was 
first introduced in [33]. Usually heaps of dimers are illustrated with the base line on the 
bottom as shown in Figure (c). If the orthogonal projection of the dimers onto the 
base line is connected the heap is connected. Further, we call dimers touching the base 
line minimal and heaps with just one minimal dimer are called pyramids. If, moreover, 
this one minimal dimer lies in the right-most column, the heap is a half-pyramid. It is 
not difficult to see that heaps of dimers corresponding to CDTs of all types of boundary 
conditions are connected. Moreover, for the type of boundary conditions considered above 
the corresponding heap of dimers is a half- pyramid, as shown in Figure ITol (c). This leads 
to a one-to-one correspondence between CDTs of topology [0, 1] x S 1 and half-pyramids of 
dimers. One can therefore enumerate CDTs and obtain the discrete finite-time propagator 
G\(x,y;t) using the generating functions for connected heaps of dimers jSHlEZj. Hence, 
the solution of the discrete problem in two-dimensional CDT reduces to a one-dimensional 
counting problem. 

There exists another interesting bijection between half-pyramids of dimers and Dyck 
paths. Consider a regular square lattice in the positive quadrant in the plane. A Dyck 
path is a one-dimensional random walk which starts at zero and finishes at zero without 
leaving the positive quadrant and only consisting of north-east steps f and south-east 
steps \ (Figure El (d)). The number of Dyck paths of length 2n is given by the n th 
Catalan number C„. 



4.3 Causal dynamical triangulations as heaps of dimers and Dyck paths 



45 




i — 

(d) 
00 

Figure 15: (a) The tiling of diamonds corresponding to the triangulation shown in Figure 
PUT! (a) is obtained by extending the triangulation to a new vertex below and above and 
removing all space-like links, (b) Starting from the dual triangulation (Figure El (b)) 
one obtains a half-pyramid of dimers by removing all space-like edges and pushing all 
remaining time-like edges (dimers) to the left until they touch each other, (c) Illustration 
of a half-pyramid of dimers with the base line on the bottom. One sees that only the 
right-most dimer is minimal, (d) Dyck path corresponding to the half-pyramid of dimers 
shown in (b). There is an one-to-one correspondence, where each up- wards step in the 
Dyck path corresponds to dimer in (b) with the respective ordering, (e) To obtain the 
original triangulation from the Dyck path one places a diamond (pair of triangles) under 
each upwards step of the Dyck path. Then, gluing together the triangles horizontally and 
removing the upper and lower row of triangles, one recovers the original triangulation. 
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5 CONCLUSION 



The correspondence between half-pyramids of dimers and Dyck paths is illustrated in 
Figure Hoi (d). Each dimer of the half-pyramid corresponds uniquely to an upwards step 
in the corresponding level of the Dyck path with the same relative ordering as in the half- 
pyramid. This gives a one-to-one correspondence between half-pyramids of dimers and 
Dyck paths. 

To obtain the original triangulation from the corresponding Dyck path one places a 
diamond (pair of triangles) under each upwards step in the Dyck path. Gluing together 
these triangles horizontally and removing the upper and lower row of triangles, one recovers 
the original triangulation. 

In summary, we have shown the following bijections 

CDT <-> half-pyramid of dimers <-> Dyck path, (139) 

which relate two-dimensional CDT to other well-known one-dimensional combinatorial 
structures. 

5 Conclusion 

In this review we have described the motivations and concepts behind the recent approach 
of Causal Dynamical Triangulation. In the case of two dimensions the model is analytically 
solvable as we demonstrated by explicit calculations. After the solution in the discrete 
setting was obtained, we showed the existence of a continuum limit, leading to a well- 
defined continuum theory of two-dimensional Lorentzian quantum gravity. We calculated 
several continuum quantities such as the effective quantum Hamiltonian, its spectrum 
and eigenfunctions, the finite time partition function and the loop-loop correlator. Using 
these results we obtained physical observables such as the expectation value of the length, 
higher moments and the effective Hausdorff dimension. In terms of these observables the 
resulting theory described a fluctuating two-dimensional quantum "universe". Further 
we have seen in this context that two-dimensional quantum gravity with Euclidean and 
Lorentzian signature are distinct theories. 

In the second part of the review we discussed the possibility of including a sum over 
topologies in the path integral. The presence of causality constraints imposed on the path- 
integral histories enabled us to derive a new class of continuum theories by taking an un- 
ambiguously defined double-scaling limit of the bare cosmological and Newton's constant. 
These causality constraints were crucial to resolve the problem of super-exponential growth 
in the number of configurations which usually appears in the context of two-dimensional 
Euclidean quantum gravity with topology changes. The resulting model of two-dimensional 
Lorentzian quantum gravity with topology changes has, besides the well-known geomet- 
rical observables, new "topological" observables, such as the finite space-time density of 
wormholes. Further, we observed that the presence of wormholes in our model leads to a 
decrease of the "effective" cosmological constant, reminiscent of the suppression mechanism 
considered by Coleman and others in the context of a Euclidean path integral formulation 
of four-dimensional quantum gravity in the continuum. 
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In the last part we discussed aspects of universality of the two-dimensional model. We 
saw that including a higher curvature term in the action does not affect the continuum the- 
ory. Further, this higher curvature term could be interpreted as be equivalent to including 
squares in the space-time triangulation. This provided evidence that the same continuum 
theory can be obtained unambiguously without dependence on the details of the discrete 
building blocks. Thus two-dimensional Lorentzian quantum gravity defined through CDT 
is a well-defined nonperturbative continuum theory. Further, we established certain one- 
to-one correspondences between CDTs, half-pyramids of dimers and Dyck paths. This 
provided a relation between the counting problem appearing in the discrete setting of 
CDT and certain well-known one- dimensional combinatorial structures. 

In conclusion, we have seen that Causal Dynamical Triangulations is a promising non- 
perturbative approach to quantum gravity. It might be the right candidate for the chal- 
lenging quest of quantizing gravity. 
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A Lorentzian angles and simplicial building blocks 

In this appendix, a brief summary of results on Lorentzian angles is presented, where we 
follow the treatment and conventions of [58J. 

Since in CDT one considers simplicial manifolds consisting of Minkowskian triangles, 
Lorentzian angles or "boosts" naturally appear in the Regge action as rotations around 
vertices. Recall from Section 12.21 that the definition of the Gaussian curvature at a vertex 
v is given by (JHJ), 

K v = (A.l) 

where e v = 2n — J2i^ v ®i ls ^ ne deficit angle at a vertex v and V v is the dual volume of the 
vertex v. Recall that the space-like deficit angle e v can be positive or negative as illustrated 
in Figure El Furthermore, if the deficit angle is time-like, as shown in Figure EE it will 
be complex. The time-like deficit angles are still additive, but contribute to the curvature 
(jA.ljl with the opposite sign. Hence, both space-like defect and time-like excess increase 
the curvature, whereas space-like excess and time-like defect decrease it. 

The complex nature of the time-like deficit angles can be seen explicitly by noting that 
the angles between two edges a, and bi (as vectors in Minkowski space) are calculated 
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A LORENTZIAN ANGLES AND SIMPLICIAL BUILDING BLOCKS 




Figure 16: Illustration of a space-like (a) and a time-like (b) Lorentzian deficit angle e v at 
a vertex v. 



using 

(a; A) y/(Si,ai)0iX) ~ (SiJi) 2 

cosOi = ; =; ^ i , sm^ = - 1 _ _ i , (A.2) 

{St, Si) 2 (bi, hi) 2 (a i; 5i) 2 (6i, 6<) 5 

where (■, •) denotes the flat Minkowskian scalar product and by definition, the square roots 
of negative arguments are positive imaginary. 

Having given a concrete meaning to Lorentzian angles, we can now use (JA.2J) to calculate 
the volume of Minkowskian triangles which we will then use to explicitly compute the 
volume terms of the Regge action. 

The triangulations we are considering consist of Minkowskian triangles with one space- 
like edge of length squared I 2 = a 2 and two time-like edges of length squared I 2 = —aa 2 with 
a > 0. The general argument a > is used to give a mathematically precise prescription of 
the Wick rotation, but it can be set to a = 1 after the Wick rotation has been performed. 
With the use of (jA.2j) we can calculate the volume of such a Minkowskian triangle, yielding 

a 2 , 

Vol (triangle) = — y/Aa + 1. (A. 3) 

Now one can define the Wick rotation W as the analytic continuation of a i— > —a 
through the lower-half plane. One then sees that for a > | under this prescription 
% Vol(triangle) i— > — Vol(triangle) (up to a 0(1) constant which can be absorbed in the 
corresponding coupling constant in the action). This ensures that 

yy . giSRcggcO^ 01 ') |_» g-^ReggeCT 6 ")^ a > (A. 4) 



In the following we set a = 1 again. Generalizations of this treatment to dimension d = 3, 4 
can be found in 0. 
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B Different boundary conditions and generating func- 
tions 

In this section we want to investigate different spatial boundary conditions of the simplicial 
manifolds considered in the construction of CDT. As we have seen in Section E3J the key 
to the solution of the model is knowing the eigenvalues of the transfer matrix, 

G x (l m ,l out ;t = l) = e- Xa2 ^ +l °^ Yl X > 

TeW(T): 

where all triangulations T are single strips of height At — 1. Considering a certain bound- 
ary condition on the simplicial manifold then reduces, in terms of the transfer matrix, to a 
simple combinatorial counting problem of evaluating the number of all possible configura- 
tions of a single strip with the same boundary conditions. The results for strips of different 
boundary conditions can be summarized in the following expression [59J: 

l^in l^out / "I _i_ 7 \ 

n x ,(l. ] l jn W / Hn T lout\ -Xa 2 (l ln +l out ) / R n\ 

^X,k{hn, I'out, t — -U — /, , , kr(] , _ 1 \k rr \ / i e ' ^ ' 

y^in \ ^out) \ h m ~T~ h out -v \ L m / 

where k denotes the multi index k = (k in , k out , k r , k rr ) e {0, l} 4 . Let us now discuss which 
cases of boundary condition we have considered: 

1. Topology [0, lJxS' 1 (as already considered in Section corresponding to k r = 1 and 
k rr = 0. Furthermore, one can mark (not mark) the initial loop by setting k in = 1 
(kin = 0) and respectively for k out for the final loop. 

2. Topology [0, l]x[0, 1] without restrictions on the orientation of the left- and rightmost 
triangle, corresponding to ki = 0, for all i G {in,out,r,rr}. 

3. Topology [0, 1] x [0, 1] without restriction of the triangle orientation on one side and 
with one up pointing triangle on the other side, corresponding to k out , k rr = and 

kin, k T 1 . 

4. Topology [0, 1] x [0, 1] without restriction of the triangle orientation on one side and 
with one down pointing triangle on the other side, corresponding to ki n , k rr = and 
kout, k r 1. 

5. Topology [0, 1] x [0, 1] with one up pointing triangle on the one side and one down 
pointing triangle on the other side ("staircase"), corresponding to A;, = 1, for all 
i G {in, out, r, rr}. 

Another useful way of characterizing cases (l)-(5) is by their generating functions. 
Hereby one can take advantage of graphical methods, where one assigns a factor of gx to 
every up pointing triangle "A" and a factor of gy to every down pointing triangle "v"- 
Considering the same cases as above: 
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B DIFFERENT BOUNDARY CONDITIONS 



1. The generating function of the propagator for triangulations of topology [0, 1] x S 1 
with one mark on the initial loop we have already obtained in Section 1231 namely 

an, r , g -, l) = lr -^l— y (B,) 

and analogously for one mark on the final loop by interchanging x and y. The 
corresponding propagator with both loops marked then reads 

[l-gx-gy) 2 

2. The generating function of the propagator for triangulations of topology [0, 1] x [0, 1] 
without restrictions on either side is in most of the literature used in its "degenerate" 
form (where initial and final boundary of length zero are allowed) , 

G opcn (x, y-g- 1) = 1 + A + V + - = V (A + V )' = z • (B.5) 

^ l-gx-gy 

3. The case of open boundary on one side and one up pointing triangle on the other 
side can then easily be generated by 

GT n (x,y;g; 1) = A • G°^(x,y;g; 1) = 9 - , (B.6) 

1 - 9 x - gy 

which includes the zero loop in the final spatial boundary. 

4. The case of open boundary on one side and one down pointing triangle on the other 
side is respectively generated by 

G°^(x, y- g- 1) = V • G°^(x, y; g; 1) = ™ , (B.7) 

v l-gx-gy 

which includes the zero loop in the initial spatial boundary. 

5. The generating function of the propagator for triangulations of topology [0, 1] x [0, 1] 
with one up pointing triangle on one side and one up pointing triangle on the other 
side is given by 

y; g; i) = a ■ G°^(x, m g; 1) • v = 1 ^ , (B.8) 

l-gx-gy 

which has obviously no degeneracies. 

Triangulations with staircase boundaries can be related to triangulations of topology 
[0, 1] x S 1 with both spatial boundaries marked [53] . 



51 




Figure 17: Construction of a triangulation of topology [0, 1] xS 1 with both spatial bound- 
aries marked by a composition of two triangulations with staircase boundary (I) and (II). 



Therefore consider two triangulations with staircase boundary, one which has the se- 
quence of up pointing triangle on the left (I) and one which has the sequence of down 
pointing triangles on the left (II) (Figure[T7J). The triangulation of topology [0, lJxS* 1 with 
both spatial boundaries marked can then be obtained by gluing triangulations (I) and (II) 
to a cylinder, where one superimposes the respective left and right most triangles. Both of 
the sequences of up and down pointing triangles uniquely define both of the markings on 
the initial and final spatial boundary (Figure IT7|) . In terms of generating functions of the 
one-step propagator, the correspondence can easily be seen by noting that the triangula- 
tions (I) and (II) both have the same one-step propagator (jB.8|) . The one-step propagator 
of the composition reads 

^in,out (* C 5 Hi 9i 1) — „2™ (^staircase ( X i Vi 9i -Q) 

9 -Ey 

= fr± , (b.9) 

(1 -gx- gy) 2 ' 

where the division by g 2 xy is due to the overlap of the left and right most triangles of 
(I) and (II). The composition (jB.9|) is sometimes denoted as two-seamed correlator. A 
generalization of this treatment, where one considers more than two strips (of either type 
(I) or (II)) in the construction of a cylindrical triangulation can be found in 



C An alternative derivation of the Hamiltonian 

In this appendix an alternative derivation of the Hamiltonian is presented, where in contrast 
to Section 12.51 the continuum limit is obtained in length space instead of the Laplace 
transformed "momentum space" . 
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C AN ALTERNATIVE DERIVATION OF THE HAMILTONIAN 



Starting point of the derivation is the one-step propagator for general boundary condi- 
tions (in the Euclidean sector), as presented in Appendix iBl 

g^.fa;« = i) = <fal*HO = ft + h) V +h _ 1)fc , ( <1 t'') <^> 

To obtain the continuum limit we use the following canonical scaling relations 

L\ — l\ • a, L 2 — l 2 ■ a, T = t ■ a, (C.2) 

where the physical correlation lengths L 1; L 2 and T stay fixed under the simultaneous limit 
a — > and 1%, I2, t — > 00. Further, one has to introduce an additive renormalization for the 
bare cosmological constant A (cf. Section l2*3|) 

A = ^ + A, (C.3) 
cr 

where A denotes the renormalized cosmological constant. The continuum propagator can 
then be obtained by 



G\ t k(Li, L 2 ; T) ^ aP G A= iog2 +A k (h = — , h — —\ t — — J 



(C.4) 



where a v is a wave function renormalization. To extract the Hamiltonian we consider the 
infinitesimal time evolution of the wave function, 

poo 

^(L 2 ,T + a) = / dv{L 1 )G Aik {L 1 ,L2-,T = a)i>{L 1 ,T) 
Jo 

POO 

= / dfi{L!) (L 2 \ (1 -aH + 0(a 2 )) \L l )^{L l ,T), (C.5) 
Jo 

with the measure dfi(Li) = V{dL\. For the further calculation it will be useful to introduce 
a new length variable D: = L 2 — Li. In terms of this variable we can rewrite (jC5j) as 



iP{L 2 , T + a)= j 2 dD(L 2 - Df G A>k {L 2 - D, L 2 ; T — a) ^(L 2 — D, T). (C.6) 

In the following we want to determine G Ajk (L 2 — D, L 2 ;T = a) up to first order in a. 
Therefore we start expanding (|C.1|) for large U and introduce new variables I := h + h 
and d := l 2 — h- We first consider the binomial factor in (jC.l|) making use of Stirling's 
approximation n\ = y/2irn (^) n (l + t|- + 0(^)1 jHO], hence 

^ . = fh + h\ e „ A a 2 (z 1+ i 2 ) = l } e -c?\i 



'1 / (l(l-f))'(l(l + f))' 

VJ 21(1 - ^ r * f 1 + 15 - - «aW + ^ d) ' (a7) 
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with 



F(/,d)=(l-j) 1 + 7) = e ^(M)io g (i-f)4(Hf)^K). (c.8) 



Introducing the scaling relations ()C.2|) into (|U.8|I . one can see that to lowest order it behaves 
like F(L 2 , D) ~exp(— ^j)- Since we want to obtain a delta function to lowest order, i.e. 
F(L 2 , D) ~ 5(-D), we can conclude that D must scale as D ~ i/a as a — > 0. Keeping this 
in mind, we can introduce the scaling relations ()C2|) - (jC3|) into (jC7|) and collect all terms 
up to first order in a, yielding 

A = x l^^e 4aL 2 ( 1 + — — 2aAL 2 - 



nL 2 V 8aL2 



a 7 10 L> 4 D 6 s . 



8L 2 32 L 2 , 96 aL\ 128a 2 L* 

In an analogous way one can also expand the remaining part of the propagator including 
the measure factor and the wave function renormalization factor, 

b ■= aPil^n^ih + hy^ih + h-iy^ 

-f) I 1 + ^(K + k rr ~ 2k in - 2/i) + 

a D 2 

4~ k rr —-z h 2 [2(^r 4" &rr 4" k r k rr ) + 4(/uj n + [l){ki n 4~ /i 1) 
2i>2 8i>2 

- 4(£v + fc rr )(fc in + /i)] + 0(a§)}, (CIO) 

where we used that ki = kf. Having obtained the continuum propagator up to first order 
in a, we can now write (10.61) in the limit a^Oas 



dDA(L 2 , D) B(L 2 , D)il—D-^- + — -^-^ + 0(a^)j ^(L 2 ,T). (C.ll) 

Upon introducing A(L 2 , D) and B(L 2 , D) as obtained above, one can then simply perform 
the gaussian integrals in (jC.llj) by doing the variable substitution D i— > /) = D/y / 2aL 2 . 
The final result reads 

-^J (l-aF + C(ai)J ^(^2,T), (C.12) 

where the Hamiltonian is given by 
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with the coefficients 

ci = 1 + 2{k in + fi) - k r - k„ (C.14) 

c 2 = (kin "I" A 4 ) {kr ^rr)(^in "I" A 4 ) "I" —k r k rr + — k rr . (C.15) 

From (jG.12j) one can further see that to obtain the right behavior to lowest order we must 
have 77 = — 1 and /i = k r + k rr — k in — k out (the remaining constant factor can be absorbed 
into the measure). For the Hamiltonian (|G.13|) to be selfadjoint on the Hilbert space 
7i = C 2 (R + , V^dL^ with fi = k r + k rr — k in — k out , we must symmetrize the Hamiltonian by 
setting k in = k out . 

Consider for example the case of topology [0, ljxS* 1 , corresponding to k rr = and k r — l; 
formula (jC.12j) then reduces to the Hamiltonian for the case with one marking on the initial 
spatial boundary, i.e. k in = l, 

H(L, — ) = —L-— + 2AL, dix(L)=L~ 1 dL (top. [0, 1] x S\ marked) (C.16) 
and the unmarked 

H(L,—) = -L-—-2— + 2AL, dfx(L) = LdL (top. [0, 1] x S 1 , unmarked), (C.17) 

which coincides with the results obtained in Section 12.51 Another Hamiltonian one en- 
counters frequently is the one belonging to topology [0, 1] x [0, 1] without restrictions on 
the time-like boundaries, i.e. fcj = 0, 

A{L 'dL ) = ~ L dtf~dL +2AL ' dfi(L) = dL (top. [0,1] x [0,1]), (C.18) 

One sees that the Hamiltonian for open boundary conditions is very similar to the one for 
cylindrical boundary conditions. Further, the spectra are the same up to some difference 
in the zero point energy (cf . Appendix El) • Therefore, one often uses the open boundary 
conditions instead of cylindrical boundary conditions for simplicity (as we partly do in 
Section EI). 



D Generalizations and the Calogero Hamiltonian 

In this appendix we want to analyze the properties and spectrum of Hamiltonians which 
one usually encounters when studying two-dimensional Lorentzian quantum gravity. As 
shown in |61j . these models can be related to the well-known Calogero model. 

In the previous appendix we have already seen that for different boundary conditions 
and different topology of the spatial slices the Hamiltonians differ slightly, but are generally 
of the form 

^4>=-^- c 4 +2Ai (d-d 
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where c\ and C2 are real constants. In the following we therefore want to investigate this 
generalized class of Hamiltonians. 

It is immediately checked that the Hamiltonian (|D.1J) is selfadjoint on the Hilbert space 
H = C 2 {R + , dfi{L)), where the measure dfi(L) is given by 17 

d/i(L) = L^dL, n = — - 1. (D.2) 

Further, for the boundary components of the partial integration in {Hip \ (/)) = {ip \ H<p) to 
vanish, the wave functions on the Hilbert space must satisfy 



(D.3) 



To have a Hilbert space with flat measure, one can pull the measure into the wave 
function by substituting ^f(L) = L^iip^L), where i/j(L) is the wave function corresponding 
to (jD.l|) . Commuting L~z through the Hamiltonian gives 

n ^ ) = - ClL ^- c ^ L+2kL+ (£!__£?)! 1, ML) = iL (D .4) 

This pulling in and out of the measure is similar to what one does when introducing 
markings on the boundary loops. 

We now consider the eigenvalue problem of the Hamiltonian (jD.4|) . 



- C ^- C 4 + 2AL+l£l 4^Z- £ )*< L > = - < D ' 5 > 
Let us perform a change of variables and wave functions, 

£ = f^ 2 , ^) = v^(y), (D.6) 

where the latter guarantees a flat measure d[i(ip) = dip for . The eigenvalue problem 
then reads 

Id 2 loo 1 A 



„ + -coW z-E] $(<p) = 0, (D.7) 

2 dtp 2 2 Y 8</? 2 J K ' 

where we set uj = y/2c\A and A = 1 — 6 (— oo, 1). This is nothing but the eigenvalue 
problem corresponding to the Hamiltonian of the one-dimensional Calogero model, 

kc ^'°^--\w + r^- l 4' M * ] = i,p - (D - 8) 

Note that the parameter range A = l— 4/i 2 G(— oo, 1) is the maximum range for which the 
Calogero Hamiltonian is selfadjoint. Two-dimensional Lorentzian quantum gravity with 

17 Note that in [STJ the definition of fi differs by one to the one used here. 
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open boundary conditions corresponds to the value A — 1, whereas for circular boundary 
conditions it corresponds to A = —3. In Appendix |E] an integrable model is introduced 
which continuously covers the parameter range — 3^A^1. The first connection between 
2d CDT and the Calogero model was established in [HT] , where a generalized CDT model 
was introduced which continuously covered the parameter range O^A^l. 

The spectrum of the Calogero Hamiltonian is well-known, further it can be related to 
the spectrum of the radial solution of the three-dimensional Schrodinger equation with 
potential U(r) = D x /r 2 + -D2?" 2 |02j- Explicitly, the solution of the eigenvalue problem 
()D.7|) can be obtained by doing a variable transformation L = uap 2 , yielding 

ii + I^ l £+ * + ^W)=0. (D.9) 



8L 2 291 4 2oo IQL 

Note that L is just a scaled version of L, so we could have directly arrived here from (JUT 
using the wave function transformation from ^(L) to &(L). From ()D.9j) one can see that 
asymptotically for L — > 00 the solution should be proportional to exp(— L/2) and for small 
L — > it behaves like L 1//4+At,/2 . Note that these are the only asymptotics compatible with 
the requirement (jD.3|) on the wave functions. Hence, we make the following ansatz for the 
wave functions $(L) oc exp(-L/2)L 1/A+fl/2 r](L). Substituting this into (|D.9J) yields 

Lr)"(L) + (1 + /* - L)rf(L) + (JL.-lfr+ 1)) nil) = (D.IO) 

This equation is known as Kummer's equation JUU], whose solutions are the confluent 
hypergeometric functions, 

V (L) = x F x {-n-l + n-L), (D.ll) 

where n—^- — |(1 + //) must be a nonnegative integer. In this case the power series of the 
confluent hypergeometric function is truncated to a polynomial of degree n, namely, the 
generalized Laguerre polynomials 



with 



. r(n + i)r(u + 1) 1 / \ , 

x F x (-n; fj,+ l;z) — ^(n + ^+i) ( D ' 12 ) 



1 d n 

I£(z) = -e z z^—(e- z z n ^) (D.13) 
" v ' n\ dz n v 



n 



^ ^ T(n - k + 1)I> + k + 1) fc! ' ^ D ' 14 ' ) 



A: 



/c=0 



yielding the final expression, 



x F x (-n; a* + 1; z) = £(-1)* (f) ... , 1W „ ,^ ... I u IT- (D.15) 



(/x + l)(/x + 2)...( A t + A;) fc!' 
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Since n— ^ — + must be a nonnegative integer, the spectrum of the Calogero 
Hamiltonian H Calo9ero reads 

E n = u(2n + fi + 1), n = 0,1,2,... (D.16) 

The corresponding eigenfunctions are given by 

$ n (<p) = C n e"^V +M iFti-n; 1 + tV), <2//(» = d<p. (D.17) 

Further, since the Laguerre polynomials are orthogonal functions, the wave functions form 
an orthonormal basis of the Hilbert space, where the normalization constant is given by 



where we used the following orthogonal relation for the generalized Laguerre polynomials 



dLe- L L»K{L)W m {L) = { ° r(M+n+1) *° rn ^ m ' (D .19) 

1 ^ ' for n = m, 



T(n+1) 

which is valid for fi ^ . 

Let us now come back to the original problem, analyzing the Hamiltonian 

The spectrum is obviously the same as for the Calogero Hamiltonian, hence from (jD.16|) 

we get 

E n = v / 2c i A(2n + /i + 1), n = 0,1,2,... (D.21) 

Further, the eigenfunctions of H(L, 8l) can be obtained from (jD.17|) by a simultaneous 
variable and wave function transformation, yielding 

_ /Mr [2A 

i/; n (L) = Ane y ^ 1 F 1 (-n;l + /i;2J—L), dfi(L) = L^dL (D.22) 
where the normalization factor is given by 



, „s -*£ r (M\" V I r(n + f i + l)~ 

A " = 2 ' Cl C " = (-J V r (" + i)r(/. + D' - ( D - a) 

Having obtained the eigenfunctions for the class of Hamiltonians H(L, di), we are now 
able to give an explicit expression for the (Euclidean) finite time propagator or loop-loop 
correlator, as defined in Section l2~5| 

G A (L 1 ,L 2] T) = (L 2 \e- Tit \L 1 ). (D.24) 
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Since the eigenfunctions form an orthonormal basis of the Hilbert space, we can insert the 
unit operator / = ^ n>0 | n) {n \ , yielding 

oo 

G K {L U L 2 -T) = ^2(L 2 | n)e~ TEn {n \ L x ) 

n=0 

oo 

= J2 e ~ TEn ^(L2)MLi) (D.25) 



n=0 



Upon inserting the energy eigenvalues (|D.21|) and eigenfunctions ()D.22)) into (jD.25)) one 
obtains 

\ Cl J ^ T{n + l)r(/i + l) 2 

/2A /2A 
x^-n; 1 + ^ 2 J — L x ) iF 1 (-n; 1 + ^ 2 W — L 2 ), (D.26) 

where we used the notation z = e~ 2v " 2ciAT . To evaluate the above summation we make use 
of the following quadratic relation satisfied by the confluent hypergeometric function 

00 T(n + fi + 1) 

^ ^ r(n + i)r(/i + D 2 1 + & ^ 1 + w v) = 



n=0 



where denotes the modified Bessel function of the /x-th kind. Gathering all terms 

together leads to the final expression for the (Euclidean) finite time propagator 

Gk{LuL 2 ,T) = \ —(LiL 2 )~2 " . w =- I„ ( 4= .T A ^ ) • (D.28) 



ci sinh(VclAT) \ V^T sinh( VciAT) 

For /i = 1, (|D.28|) corresponds to the "unmarked" propagator for the 2d CDT model with 
circular boundary conditions, as obtained in (JoTj). where the identity was given by 

1 f°° 
(Lt I L 2 ) = —8{L X -L 2 ), / dL I L)L(L\= 1 (D.29) 

L i Jo 

Note that by changing the measure one can always shift around factors of L's from f)D.28|) 
to (JD~29|) . 

In Section 12.61 we were interested in calculating expectation values of the length and 
higher moments with respect to a certain eigenstate \ n), 

POO 

(L m ) n = (n\L m \n)= / d^{L)^ n {L) L m ^ n {L). (D.30) 
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The result can also be obtained for the generalized class of Hamiltonians IjD.lj) . as consid- 
ered above. Inserting the eigenfunctions (|D.22|) into ()D.30|) yields 



(2/*) = (i±^ f 1 ( " ' 1 



dL L mH 



AAJ r(n + /x + l) 
For the evaluation of the integral we use the following integral expression [U4*j . 



(D.si; 



dLL^ +m e-' L 



[1 + fi) n (-m) n + m + 1) 



x 



r(n + l) 2 

x 3_F 2 (— n, // + m + 1, m + 1; 1 + //, 1 + m — n; 1), // > — m — 1, (D.32) 
where 3 F 2 (ai, a 2 , 03; 61, b 2 ; z) are the generalized hypergeometric functions, defined by 



pFg(ai, a p ; 61, b q ; z) ^1 pjP (& ) fc ' 

and (a) n are the usual Pochhammer polynomials, 

. T(a + n) 
r(a) 

Hence, inserting ()D.32|) into ()D.31|) gives the final expression for the moments, 



(D.33) 



(D.34) 



(L m ) = — 

\ /n V A A 



ci \ f r(n — m)T(/i + m + 1) 



x 



AAJ r(n + l)r(/i + i)r(-m) 
x 3 F 2 (-n, + m + 1, m + 1; 1 + /i, 1 + m - n; 1). 



(D.35) 



Note that the poles of T(—m) cancel with those of the hypergeometric function leading to 
a finite expression. Further, it is easy to see that all moments scale as (L m ) n ~A _: ¥. 



E Discarding unphysical double scaling limits 

In this appendix we discuss certain double scaling limits of the CDT model with topology 
change, which we discarded as unphysical in Section 13.41 In terms of the parameters a 
and f3 in the scaling relation (|95j) . these were the scalings with (3 = 1, instead of the scaling 
with (5 = I which led to a physical sensible continuum theory, as discussed in Section 13.41 
Hence, for (3 = 1 the scaling relation (J95|) reads 

h = ^h ren aVA a {X + Y) 1 - , (E.l) 

where the normalization factor on the right-hand side has been chosen for later convenience. 
To obtain the effective quantum Hamiltonian we follow the same procedure as the one used 



60 



E DISCARDING UNPHYSICAL DOUBLE SCALING LIMITS 



in Section El 41 Upon inserting (|9() )1 -()92 |1 and (jE.lj) into the time evolution of the discrete 
wave function (|96|). and expanding both sides to first order in a one obtains 

/ io ° d7 
— {A(X, Z) + B(X, Z)a + 0{a 2 )}i{j{Z), (E.2) 
-ioo 27TI 

where the leading-order contribution is given by 



Aix ' z > = { z-x)(i + c { x,z)) (E ' 3 > 



with 



C(X, Z) = y/1- h 2 ren (X - Z)~ 2 «A«. (E.4) 

For the Laplace transform of A(X, Z) to yield a delta-function, the scaling should be chosen 
such that a ^0. Consider now the terms on the right hand side of (jE.2|) . which are of first 
order in a, 



B(X Z) = K en (X + Z - AZ^A* 

1 ' ' (x - zy+ 2a c(x,z){i + c(x,z)) 2 

XZ -2A + 7 (X - Zf 



(E.5) 



(X - Z) 2 C(X, Z) (1 + C(X, Z)) ' 

one finds that for a ^ — 1 the continuum limit is independent of any new coupling associated 
to the contribution of the wormholes and therefore leads to the usual Lorentzian model. 
This becomes clear when one expands the last term of (jE.5|) in (X — Z), which yields 

XZ 2A ' XZ 2A f l + I hl n A«(X - Z)- 2 « + 0((X - Z)-^)\ . (E.6) 



(X-Z) 2 C(1 + C) 2(X-Z) 2 \ 4 

For a ^ — 1 the term depending on h ren does not have a pole and therefore does not 
contribute to the Hamiltonian. Since we are only considering non-fractional poles, this 
leaves as possible a- values only a = and a = —\. 



E.l The case f3 = 1, a = 

For a = the Hamiltonian retains a 7-dependence contained in the first line of ()E.5j) . 
Since there is no immediate physical interpretation of 7 in our model, it seems natural to 
choose 7 = 0, although strictly speaking this does not resolve the problem of explaining 
the 7-dependence of the continuum limit. Setting this question aside, one may simply 
look at the resulting model as an interesting integrable model in its own right. 



In order to obtain a delta-function to leading order, one still needs to normalize the 
transfer matrix by a constant factor 2/(1 + s), with s := ^1 — h 2 en . After setting 7 = 
and performing an inverse Laplace transformation, the Hamiltonian reads 

d 1 / d 2 d \ 



E.2 The case (3 = 1, a = — \ 61 



2 

It is self-adjoint on the Hilbert space 7i = C 2 (Jt + ,L s ~ 1 dL). Further setting L = one 
encounters the one-dimensional Calogero Hamiltonian (cf. Appendix IDj) 

rr, 9 n I d 2 1 2 2 1 A _ o , 

^•^ = -2fl? + 2^-87' (K8) 

with u; = and A = l — 4(1 — s) 2 , which implies that the model covers the parameter 
range —3 ^ 1. The maximal range for which the Calogero Hamiltonian is self-adjoint 
is —oo<A^ 1. It is interesting to see that the usual Lorentzian model without topology 
changes corresponds to A = 1 in the case of open boundary conditions, and to the value 
A — —3 for circular boundary conditions (cf. Appendix [D]). This selects the parameter 
range —3 ^ A ^ 1 as a natural choice. The Hamiltonian (|E.8|) has already appeared 
in a causal dynamically triangulated model where the two-dimensional geometries were 
decorated with a certain type of "outgrowth" or small "baby universes" [HI] . This model 
covered the parameter range ^ A ^ 1. 



The spectrum of the Hamiltonian (jE.7|) is given by 



V2A 

E n =- (2n + s), n = 0,1,2,.... (E.9) 

s 

The corresponding eigenfunctions read 

ip n (L) = Ane-^ 1 iF^-n, s, 2V2AL), d/i(L) = U~ x dL, (E.10) 

where \F\{— n, a, b) is the confluent hypergeometric function, as defined in (|55jl. The 
eigenfunctions {tp n (L),n = 0,1,2,...} form an orthonormal basis with the normalization 
factors 

A " - < 8A > V r(« + i)-rU 

Let us for completeness also state the finite time partition function or loop-loop correlator 



with s G [0, 1]. One sees explicitly that the case s = 1 or, equivalently, A = 1 corresponds 
to the pure two-dimensional CDT model with open boundary conditions. 

E.2 The case f3 = 1, a = -\ 

For a = — \ the result does not depend on 7 and therefore on the detailed manner in which 
we approach the critical point. However, the Hamiltonian 
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cannot be made self-adjoint with respect to any measure d[i{L) because the boundary part 
of the partial integration always gives a nonvanishing contribution. We therefore discard 
this possibility. 
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